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Matrix Inversion by a Monte Carlo Method’ 


The following unusual method of inverting a class of matrices was de- 
vised by J. von NEUMANN and S. M. Utaws. Since it appears not to be in 
print, an exposition may be of interest to readers of MTAC. The method is 
remarkable in that it can be used to invert a class of m-th order matrices 
(see final paragraph) with only m? arithmetic operations in addition to the 
scanning and discriminating required to play the solitaire game. The method 
therefore appears best suited to a human computer with a table of random 
digits and no calculating machine. Moreover, the method lends itself fairly 
well to obtaining a single element of the inverse matrix without determining 
the rest of the matrix. The term “Monte Carlo” refers to mathematical 
sampling procedures used to approximate a theoretical distribution [see 
MTAC, v. 3, p. 546]. 

Let B be a matrix of order m whose inverse is desired, and let A = I — B, 
where I is the unit matrix. For any matrix M, let \,(M) denote the r-th 
proper value of M, and let M;; denote the element of M in the i-th row and 
j-th column. The present method presupposes that 


(1) max|1 — A,(B)| = max]A,(A)| < 1. 


When (1) holds, it is known that 


(2) (B>)y = (I - APs = £ (A%e. 
k=O 


The Monte Carlo method to compute (B-"),; is to play a solitaire game G;; 
whose expected payment is (B-),;. According to a result of KoLMOGOROFF? 
on the strong law of large numbers, if one plays such a game repeatedly, 
the average payment for N successive plays will converge to (B-),; as 
N — o, for almost all sequences of plays. The rules of the game will be 
expressed in terms of balls in urns, but a computer would probably use a 
table of random digits. 

For 1 =i,j =n pick probabilities p;; = 0 and corresponding “‘value 
factors’’ v,;, subject only to the conditions that p,;; = ai;and "je: pis < 1. 
Let the “stop probabilities” p; be defined by the relations p; = 1 — Dja1 pi. 
Take m urns. In the i-th urn U,(¢ = 1,2, ---,m) put an assortment of 
n + 1 different types of balls. Each ball of the j-th type is marked “‘j,” 
and will be drawn from U; with probability ;;(j7 = 1,2,---,m). The 
(m + 1)-th type of ball is marked “STOP,” and will be drawn from U; with 
probability p,. 

For the game G,; now to be defined, the value of a play is a random 
variable G;; whose expectation will be proved to be (B-'),;. First draw a 
ball from U; (all drawings are with replacement). If it is a STOP ball, the 
payment G,; is 0 (if ¢ = 7) or p;~ (if i = 7). Otherwise the ball must carry 
a mark 4,(1 = 4; = m), and one is next to draw a ball from U;,, which in 
turn tells one whether to stop or draw again. One follows this treasure 
hunt from urn to urn until a STOP ball is first drawn—say from U, on the 
k-th drawing. If kh ~ j, the payment G,; is 0. If 4 = 7, suppose one has 
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arrived at U; via a route t = tp —> 41 —> i2 — «++ — %_1 Jj, whose abbrevi- 
ation is p. Then the payment G;,; is defined to be 
Vipir = Vig Ving *** Vi, gidi 

THEOREM 1. If (1) holds, Gi; has the expectation (B-),;. 

Proof. The probability of following the route p and then drawing a 
STOP ball is 

Pi = PigiPigig *** Pig sidi- 

The expected payment in the game G;,; is 


E(Gis) a os 2 (P,bi)( V bi) = ) P,V,, 


where the sum is extended over all routes p from 7 to 7. Since p;vi; = aij, 


EGis) = 85 + DO Lees Le Gading °° Cys 
k=1 i1=1 %—-1=1 


= Ty + © (AM is 


k=1 


where 6;; is the KRONECKER symbol. By (2) the theorem is proved. 

Let us calculate the variance o7, of the payment G,; about its expected 
value. Let R be the matrix (a,,0;;) = (piv7,). If max;|d.(R)| = 1, the proof 
below can be modified* bs show that oi; = ©. Assume max;,|A,(R)| < 1, 
and define (J — R)“ = 

THEOREM 2. If pt-F (R)| < 1, 


oi; = Type — (BY; 
Proof. Using the above notation 
oi, = E{Gy — (B).s}? = E(Gi,) — (Bi, 
= ¥ (P,6;)(V.p7)? — (BY) 
a 


= p77 P,V, — (B)is. 


If max;|A;(R)| < 1, the end of the proof of Theorem 1 can be modified to 
show that the last sum is 7;;. This proves Theorem 2. 

One can of course compute one whole row of B-' at once by playing the 
games G;.(a = 1, 2, ---,) simultaneously. It may even be practicable to 
play all the games {G,;} (7, 7 = 1,2, ---, m) simultaneously, using not only 
the full route 9 > 4: > i, > --- —-%1— j, but also each partial route 
dy > Sy41 > +++ —>%-1— j as a separate play. In any case the arithmetic 
can be reduced by accumulating the totals of V, for a series of plays and 
dividing by Np; at the end, to obtain the average payment after N plays. 

If a3, = 0 and >521 a; < 1, one can take p;; = a; and v;; = 1. Then 
V, = 1, and the computer has only to make a series of plays and keep a 
tally of the frequency with which the STOP ball was drawn in U; (j = 1, 
2, +++, ).4 If N is selected to be a power of 10, the only arithmetic required 
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to invert B is nm divisions of a vector by p;. When »,, = 1, A = R and 
T = B—, and hence of, = p;-(B) [1 — p;(B™)«], corresponding to the 
variance of the binomial distribution. 

GeEorGE E. ForsyTHE 


RicHarD A. LEIBLER 
NBS Institute for Numerical Analysis 
Univ. of California, Los Angeles 
2901, 18th St. 
Washington 9, D. C. 

1 The preparation of this paper was sponsored (in part) by the Office of Naval Research. 

2A. Ko_mocororF, Grundbegriffe der Wahrscheinlichkeitsrechnung, New York, 1946, 
p. 59. The writers are indebted to T. E. Harris for this reference. 

3 The fact that o;; = © does not interfere with the convergence of the average value 
of N games to (B~),;;. However, conventional error estimates in terms of variances no 
longer apply and, in at least certain matrix inversions where oj; = ©, the accumulated 
payment after N games cannot be normed so as to be asymptotically normally distributed 
as N—> o. See W. Fetter, “Uber den zentralen Grenzwertsatz der Wahrscheinlichkeits- 
rechnung,”” Mathematische Zeitschrift, v. 40, 1935, p. 521-559 and v. 42, 1937, p. 301-312, 
ond "Uses = Gesetz der grossen Zahlen,” Szeged, Acta Univ., Acta Scient. Math., v. 8, 
1937, p. 191-201. 

‘ fi is this case which we learned from von Neumann and Ulam. 


Maximum-Interval Tables 


Both the article by HErGEt & CLEMENCE [MTAC, v. 1, p. 173-176] and 
the note by MILLER [MTAC, v. 1, p. 334] on optimum-interval tables ignore 
the possibility of a continuously variable interval. It is of some interest to 
examine the reduction in the number of entries made possible by what might, 
by analogy, be termed ‘‘maximum-interval”’ tables. Using the principles of 
optimum-interval tables, with the modifications suggested by Miller, the 
tabular values of the argument are no longer restricted to terminating deci- 
mals so that the interval may be allowed to assume at each point the 
maximum value consistent with the stated allowable error. 

The chief objection to a punched-card table in this form is that all (or 
nearly all) of the digits in the argument will have to be used in the inter- 
polating factor. In some cases this objection could be overcome by inserting 
additional cards corresponding to values of the argument terminating in the 
appropriate number of zeros, or by splitting the whole range into a number 
of sub-ranges, in each of which the allowable error is varied slightly to make 
tabular arguments coincide with the end-points; or, of course, by an addi- 
tional operation of subtraction. Generally, however, the saving in cards is 
not worth the additional cost of preparation and the resulting complication. 

Herget in The Computation of Orbits [see MTAC, v. 3, p. 418-9] 
gives an optimum-interval table of x~*/? using quadratic interpolation, with 
a note “This is the first time such a table has ever been printed for use with 
a hand calculating machine.” A human computer can easily exercise the 
requisite judgment to use the continuously variable intervals of a maximum- 
interval table; in the simplest case it will involve nothing more serious than 
allowing the interpolating factor occasionally to exceed unity in a particular 
digit. A punched-card machine can only do this at the expense of a separate 
operation. There may, therefore, be a use for variable interval tables in 
computation by desk machines. 
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The chief field of use is likely to be in connection with automatic digital 
computers, since the number of operations and the number of digits in the 
multipliers is often less critical than the amount of high-speed storage 
required. With cubic and quartic interpolation the number of entries re- 
quired is very small indeed. 

Before considering the theory of these tables, it is desirable carefully to 
examine the significance of “‘error’’ as applied to values obtained by inter- 
polation. Corresponding to the function tabulated and the conditions of its 
use the error may be assessed in four ways: as an absolute or relative error 
in the tabulated function, or as relative to an absolute or relative error in 
the argument. Examples of these four cases are: 

(a) sin x, to a stated number of decimals; here the error is absolute. 

(b) x”, with a stated percentage error as is always required when multi- 
plying factors are used to extend tables for a limited range of the argument; 
here the error is relative to the tabulated function. 

(c) sin x, with an accuracy corresponding to a stated absolute error (e.g., 
0°.0001 or 1’’) in the argument. 

(d) x, with an accuracy corresponding to a stated relative error (say 1 
in 10’) in the argument. 

Each of these cases may occur and it is essential that the appropriate 
one.should be adopted if the utmost economy in number of entries is to be 
achieved. Errors arise from several sources, but all except those due to 
approximations made in the interpolating formula can be substantially and 
effectively reduced by the retention of an additional, or guarding, figure; 
it will be assumed that this is always done in the type of table to be dis- 
cussed, so that the theoretical errors of approximation may be directly 
equated with the allowable, or stated, errors. 

The TAYLor series for a function y of x at x = a + ¢ may be replaced, in 
the range — $h < t < + $h, by aseries in terms of CHEBYSHEV polynomials:! 


(1) ya + t) = a + arCi(4t/h) + a2C2(4t/h) +---+ apC,(4t/h) +--- 


where C,(~) = 2 cos (p arccos $£) is the Chebyshev polynomial of degree p 
and where 


ae si (th)? (th) 


let)” + 2@+DG+2)" 


in which Fy is the p** derivative of y at x = a. 

It is known that the truncation of (1) after the m** term will provide the 
most efficient polynomial approximation of degree (m — 1) to the function 
y(a + #) in the range — $h < t < + $h. The leading term of the maximum 
error of the approximation is clearly 2a,, which may for the present purpose 
be simplified to: 





hy nae 


(p) 4 pt2) 4... 


hn 
(3) igi y™. 


This may now be equated to the allowable error to give the maximum 
permissible value of h as a function of x. If €:, €2, es and ¢, are respectively 
the maximum absolute and relative errors allowable in y and x, then (3) 
must be equated with: 


(a). (b) egy (Cc) ey’ = (dd) cary’ 
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where the signs of the e’s are chosen so that h* is positive. 

Let z be an auxiliary variable which takes integer values corresponding 
to the values of x at which y is to be tabulated. Then, to a very good approxi- 
mation, h may be identified with dx/dz so that we obtain the following ex- 
pression for z in terms of x: 


(4) s< f * Ielde. 


This will give z as a function of x, and in due course x, y and its deriva- 
tives as functions of z from which the table can be prepared. 
For a simple power of x, 


y = x™ (m not a positive integer); 


substitution in and integration of (4) with the condition that z = 0 at 
x = 1, gives: 

(a) 2 = e-“*(n/m)b(m, n)(1 — x™/*) 
(5) (b) 2 = es-/"b(m, n) Inx 

(c) 2% = (me3)—/*nb(m, n)(1 — x'*) 

(d) 2 = (me)~"*b(m, n) In x 
where 

+ b(m, n) = [m(m — 1) --- (m — n + 1)2!**/n!}!*, 


the sign being taken to make z positive. Cases (b) and (d) are, of course, 
identical in form. These expressions enable the number of entries for a 
particular range of x to be determined at once. 

For example, in the seven-decimal table of reciprocals considered by 
Herget & Clemence the value of ¢; is taken as 5-10-’. A more realistic ar- 
rangement would arise from specifying maximum relative errors; a com- 
parable table should therefore be based on eg = ¢, = 5-10-'*/, the geometric 
mean of the relative errors at beginning and end. Case (c), of a stated 
maximum absolute error in x, may well arise in practical computation; ¢€; 
can be taken as 5-10~’. The corresponding expressions for z are: 


(a) 2 = 3n(4-10%)**(1 — x-/*) 
(b), (d) zs = 3(4-10")/* In x 
(c) s = 4n(4-10%)"/*(x"/" — 1) 


The number of entries required in a table to cover the range x = 1 to 10, 
with » = 2, 3, 4 (corresponding to linear, quadratic and cubic interpola- 
tion), are: 


n=2 n= 3 % = 

(a) 685 65 21 
(b), (d) 648 64 21 
(c) 2163 139 36 


The reduction to 685, for case (a), as compared with the 1368 of Herget 
& Clemence and the 924 of Miller is substantial. There is here little difference 
between the requirements of cases (a), (b), and (d), though (b) is certainly 
the more realistic table. The very considerable reduction effected by the 
inclusion of a quadratic term would suggest that it might well be worth 
while introducing the additional step into the interpolation. 
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As a second example, consider a table of x~*/* with a maximum relative 
error of 10~’? which corresponds approximately to that of Herget’s table. 
The number of entries (case (b)) required are: 


2 (linear) 3526. 
3 (quadratic) 204, as compared with Herget’s 316. 
4 (cubic) 50. 


n 
n 
n 


The coefficients of the interpolating polynomial of degree (m — 1) are 
obtained by expanding the individual terms of the truncated series (1) and 
rearranging as a power series in ¢. This series is then appropriate to the range 
— $h < t < $h. In order to facilitate computation it is desirable still further 
to transform the series as a power series in s where 


s = (a+ 2) — xo 


and x9 is a value of x terminating with several zero digits. It is only necessary 
to retain sufficient figures in the coefficients of powers of s to cover the range 
appropriate to ¢, provided they are rigorously consistent with the coefficients 
in the ¢-series; they should be calculated one at a time starting with that 
of s*-, which is of course the same as that of #*—. 

It will be noted that for y = x”, y in case (a) and x in case (c) are both 
polynomials in z of degree m — 1. In cases (b) and (d) both x and y are of 
exponential form; since the coefficient of In x in (5) does not change rapidly 
for small changes in m, it would be possible to use the same series of values 
of x (or of y) for a number of fractional powers. 

In his note Miller referred to the previously accepted convention of not 
modifying function values to reduce the error due to neglect of second 
differences. Actually such modification (for a slightly different reason) is in 
use in both the British and American surface and air almanacs; but in many 
cases there are sound reasons why the tabular values should be more accurate 
than the interpolates. It is possible to reduce considerably the error due to 
neglect of second differences without modifying the function values: in 
principle, the first difference is modified so that the error at the end of the 
range is numerically equal, with opposite sign, to the maximum error in the 
range. Specifically (3 — 2V2) = 0.043 of the double second difference is 
subtracted from the first difference; the error is then —0.043 of the double 
second difference at the end of the interval and +0.043 at a point 0.414 
along the interval. This error compares with 0.0625 without modification 
and 0.03125 with modification of the function. The device has obvious dis- 
advantages, but it has some valuable applications—mainly in navigational 
tables. 

D. H. SADLER 
Royal Greenwich Observatory 
Herstmonceux Castle, Sussex 

1 This development is taken, with only tight change of notation, from Miller's paper 

“Two numerical applications of Chebyshev polynomials,” R. Soc. Edinburgh Proc. Sec. A, 


v. 62, 1946, p. 204-210. It is, of course, simple to obtain the optimum expansions for small 
values of » by direct methods. 
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Planning and Error Considerations for the 
Numerical Solution of a System of 
Differential Equations on a 
Sequence Calculator 


1. Introduction.—This paper deals with the solution of a specific system 
of fourteen ordinary differential equations, 


(1) ai) = fit, ---, 21, 2), where 4=1,2,---,14. 


The system was solved for the United States Navy Special Devices Center 
on the IBM Selective Sequence Electronic Calculator. A discussion of 
basic theoretical planning and error evaluation will be presented herein. The 
complete solution process will be described in a joint report with Pau 
Brock of the Reeves Instrument Company, which it is hoped will be avail- 
able from the Special Devices Center. The present theoretical discussion 
outlines reasons for the choice of method, the detailed step-by-step process 
decided upon, the considerations which led to choice of interval length, and 
those which are concerned with the choice of scale. 

There were three types of error in the solution process, namely, the 
truncation error, the rounding error, and the error due to the presence of 
nonanalytic functions in the f’s. It is possible to obtain the immediate error 
of each type and to develop a uniform method of evaluating the complete 
effect of each such error. In the case of the truncation error, there is a reso- 
nance effect, a build-up of error, which is different in nature from instability, 
as this system is stable in the usual sense. The rounding error is treated by 
probability methods. However, the present treatment of this error is some- 
what simplified since the effect is small. The nonanalyticities require a special 
technique, which will be discussed later. The methods developed in this 
paper, with this specific example in mind, differ considerably from those 
of RADEMACHER.! 

2. Equations.—The f;’s in this system can be calculated from a given set 
of values, 21, ---, 14, and ¢ by means of the elementary operations of arith- 
metic and the functions: square root, absolute value, sine, cosine, and 
arcsine. A series expansion was used for the arcsine; Newton’s Method was 
used for the square root; and a table was used for the sine and cosine evalua- 
tion. The nonanalyticities appeared because of expressions which occurred 
in the form, |z; — k|, when certain of the dependent variables, z, pass 
through specific values. The initial values were in a region of analyticity, but 
nonanalyticities occurred eleven times during the solution process. 

3. General Method.—Consider f; as a function of time. Suppose f; has 
been evaluated at ¢ — h, t — 2h, t — 3h, and t — 4h. A cubic in time is 
passed through these values and the polynomial integrated between ¢ — h 
and ¢. The result is a value for Az;, the increment in z; between ¢ — h'and t. 
In the usual Milne method, the value, z;, obtained at ¢ is substituted{in the 
f's and a new cubic formed which goes through the values of f at t, ¢ — h, 
t — 2h, and ¢ — 3h. This in turn is integrated between ¢ — h and ¢ to yield 
a new estimate for Az;. The new value, z;, can be used in the same way to 
obtain another value of Az; and the process continued until there is no 
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significant difference between successive values of z;. A modification of this 
method is used. The first step in the above process is referred to as the 
“open”’ integration; the remaining steps are called ‘‘closed’’ integration. 
The value of Az; obtained from the open integration is a simple linear 
combination of the values of f; for t — h, t — 2h, t — 3h, and t — 4h, ie., 


(2) Ass = 32 [S5ft — h) — S9f,{t — 2K) + 37f,{t — 3h) — 9f,(t — 48) 
For a closed step 
(3) sy = [9f(t) + 19f(t — b) — Sf{t — 2K) + F(t — 3h)}. 

4. Step-by-Step Error.—For simplicity, let ¢ = 0 in the above equation. 


The usual remainder argument shows that, if P is a polynomial of the third 
degree which equals f;att = —h, —2h, —3h, and —4h, then for —h =t =0 


(4) fi — P= ferry) EE DE + E+ WE +O) 





for some x’ between —4h and 0. Suppose f;“” is a constant and integrate 
between —h and 0; then, the error, ¢;, in Az; of the open integration is 


(5) ef = .349f,0PH5, 


If the correct value for f; at ¢ = 0 were known, a similar argument for 
the closed integration step would yield the error in the closed steps, i.e., 


(6) ef = — .0264f,2PH'. 


However, in the first closed integration step for Az,;, the following expres- 
sion holds: 


(1) baer = 5 [9fi() + 19flt — B) — Sflt — Dh) + fat — 3], 


where f;(#) is obtained by substituting the values of z involving the open 
integration error. Consequently, the first closed integration step involves 
an error 
” 4 Of; 

(8) t= €f + ea 
Now, if & is sufficiently small, the second term is less important than the 
first. Consequently, if f;*” is essentially constant, ¢,° and e,“' differ in sign, 
and the values for the open integration and the first closed integration for 
Az; span the true value. The difference between these values forms an upper 
bound and a good estimate for the truncation error. 

It may be noted that the same argument shows that the error, e,, in 
the k* step can be expressed in terms of the error, ¢,“*-, in the previous 
step, i.e., 

of; 


(9) ect = ef + +r eee 
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This forms an easy basis for convergence considerations. The final con- 
vergence has an error which satisfies the equation, 


_ 9h 5. af 
24 ~ OZ 


5. Actual Method Used.—The f’s are complicated functions from the 
computational point of view, and their evaluation consumes much machine 
time. With this in mind, consider the relatively simple method given above. 
In view of equation (10) and the machine time required to evaluate the f's, 
there seems to be no point in carrying through the iterative process to the 
bitter end. In fact equations (8) and (10) suggest that considerations can be 
confined to an open and closed step. However, the price of an open and 
closed step in machine time is the same as the price of two open steps. It is 
possible to compare the error which follows from an open and closed step 


(10) ef 


€.°° = €/*. 


, ; . The latter 
can be obtained by dividing .349 in equation (5) by 16 and comparing the 
result, .022, with .0264 of equation (6). Thus, the open integration pro- 
cedure, which is simplest from the coding point of view, is the best. (This 
argument is based on the use of cubic approximations.) Even if this situation 
were reversed, it would be better to choose the open integration method since 
a smaller step interval is used, and, hence, error estimates are more depend- 
able. Similar objections may be raised to any other system which uses a 
larger h. . 

However, certain other considerations must be made. The choice of the 
interval, h, is necessarily based on the initial information given in the 
problem. Later in the run, 4 may be too large or too small. This can only be 
established by knowing the truncation error, which can be indicated by the 
difference between a closed and open step, as has been shown above. The 
procedure is to use an open integration at each step but, at every tenth step, 
to follow it by a closed integration. At each tenth step, the difference between 
the open and closed steps can be obtained from the printer in the sequence 
calculator, and one can infer from this whether the truncation error is 
tolerable. 

If the truncation error is small enough, one is justified in increasing h 
by an integral multiple. This is readily done using the previously computed 
values of f. The coding for the problem was set up in such a way that one 
could start by feeding in cards with the properly spaced values of f for four 
points. To increase h then by an integral multiple, it was only necessary to 
take certain cards previously punched by the machine and feed them in. 
It was not necessary to decrease h, but if necessary, intermediate values of 
f could have been obtained by interpolation involving four differences. 

6. Choice of h.—From physical reasoning and also from an actual con- 
sideration of the variational equation at the initial part of the run, it was 
determined that the system of equations is stable in the absolute sense. 
When an error is made, the effect of this error does not increase indefinitely, 
but, instead, dies down at a very slow rate. 

Initially, the possibility of a run to ¢ = 100 had to be planned for. Calcu- 
lations based on the initial situation and these considerations indicated that 


with interval, hk, and the error of two open steps with interval 
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the interval, h = .04, is too large and that 4 = .02 would be about right. 
But the cost of the solution for k = .02 was prohibitive. Furthermore, it 
was felt that the assumption that all the errors would build up in the same 
direction was too conservative, and, hence, it was decided to use the interval, 
h = .04. Due, however, to the resonance effect which will be discussed later, 
the errors do effectively build up in the same direction so the stated reason 
is not sound. On the other hand, f“Y does decrease as evidenced by the 
truncation error, and, while the total accuracy was not quite that planned, 
the result was adequate. 

To begin the process, it is necessary to compute the values of f’s at four 
points. These initial values were computed by finding the corresponding 
values of z by means of Taylor’s series and computing the f’s. Since it is 
relatively easy to increase h, an initial value of kh = .001 was chosen, and 
the f’s were computed for ¢ = .000, ¢ = .001, ¢ = .002, and ¢ = .003. Using 
h = .001, the computation process was carried to ¢ = .100. This resulted in 
values at ¢ = .00, ¢ = .02, ¢ = .04, and ¢ = .06, which were used to start 
the run with 4 = .02. At ¢ = .08 the values for z using the short interval, 
h = .001, and also an open step value for 4 = .02 and a closed step value 
were then available. Of course, the smaller h value yields a much more 
accurate value for z. It is interesting to observe that for every z, the accurate 
result lay between the open and closed values. 

The hk = .02 run was continued until about ¢ = 2.18, and the h = .04 
run began with ¢ = 2.10. Again an overlap was permitted for checking 
purposes in order to be sure that the change of interval was correctly made. 
It was found that the shorter interval value again lay between the open and 
closed value for the longer interval. 

The interval, h = .04, was continued until ¢ = 68, where one of the 
variables, z, ran off scale, and the solution was no longer of interest. During 
an earlier false run, an attempt was made to increase the interval toh = .12, 
but this yielded large truncation errors. Therefore, the interval, h = .04, was 
used again. During the final true run, no such increase was attempted. 

After the h = .04 run, the h = .02 run was repeated and continued in 
order to check the error of the kh = .04 run. Unfortunately, only the run 
from ¢ = 2.10 to ¢ = 2.90 was duplicated at an interval, h = .02. However, 
this provided a valuable method of checking the error theory developed later. 

7. Scaling Questions and Watches.—The registers in which a number 
can be entered contain 9 digit-positions. The scale for a quantity must be 
chosen so that the effect of truncation error can be clearly seen. On the other 
hand, some leeway against spill-over must be available. Consequently, the 
size of each quantity must be estimated in advance. 

For this estimation, a certain amount of physical information is avail- 
able. A set of assumptions, based on this information, were made concerning 
the sizes of the quantities involved, and values for all other quantities were 
estimated. However, the assumptions were not based on completely reliable 
information, and it was necessary to verify that, during the course of the 
computation, these assumptions remained valid. This was done by a system 
of “‘watches.”’ Certain of the quantities were printed by the machine, and 
others could be obtained readily from printed results. During the complete 
calculation, up to the time the calculation was stopped, the assumptions held. 
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8. Hand Computations.—The Taylor’s series for the z’s had to be ob- 
tained by hand computation. This was done in triplicate, as it was believed 
that the error rate was approximately about 1 in a 100 calculation steps. 
Since the amount of work was of the order of magnitude of 10,000 steps, 
a calculation in duplicate would have too high a probability for a duplicate 
error. 

In addition, it was also necessary to provide a hand computation of two 
steps of the machine computation for checking purposes. This was done for 
the ¢ = .004 and ¢ = .005 steps, but unfortunately many values appeared 
as zeros necessitating the provision of further material for checking pur- 
poses. In addition, a complete machine step consisting of an open and closed 
step at the later value, ¢ = 6.44, was duplicated by hand and the values 
checked against the machine values. In these checking computations, the 
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Illustration 1. Solid curve I is true solution. The broken curve is z as defined 
by equation (11). The jumps are yz. 


same rounding procedures as those of the machine were used. Hence, a 
check demands precise digital agreement to every digit which appears. 

Nevertheless, this numerical checking was inadequate to catch one error. 
It should be pointed out that a specific numerical calculation is logically of 
too low a type to constitute a complete coding check. A complete check 
should start with the coding and reconstruct the original mathematical 
relations. Thus, it is clear that a correct coding system is required. 

9. Definition of the Error.—It has been mentioned that there are three 
types of errors inherent in the solution. Associated with every error there are 
two problems. It is necessary to estimate the error at the point where it is 
made. It is also necessary to evaluate the effect of the error on those portions 
of the solution which follow in time the point at which an error is made. Let 
the vector zs = (21, «++, 214) denote the computed solution. Let (cz, ---, cts) 
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denote the correction vector which, when added to 2z, will yield the true 
solution. The computed solution exists only at the points at which a compu- 
tation was made. Let z be defined between such points by the condition that 
it satisfies the equations 


(11) Zu’ = ful, °°, 214, t) 


and is continuous at the left-hand end point of each such interval. At the 
right-hand end point of any interval, z has a discontinuity corresponding to 
the error made at the right-hand point. (See Illustration 1.) 

At each solution point the total error is made up of two parts. The first 
part, termed the rounding error, is the difference between the computed 








Illustration 2. 
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solution and the computationally perfect result of the same computation 
made with no rounding errors. The second part is the difference between the 
computationally perfect result and the result of a correct integration between 
the previous step and the present one. If the points involved in the integra- 
tion process do not span a point of nonanalyticity for the f being integrated, 
the second part will be termed the truncation error. If, on the other hand, 
there is a point of nonanalyticity, for the f being integrated, present in the 
span covered by the integration process, this error will be called a non- 
analyticity error. The distinction between truncation errors and non- 
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analyticity errors is set up purely with the situation of the present practical 
problem in mind and appears to be the most convenient in this case. 

The total correction vector at the a computational point corresponding 
to the total error made at that point is denoted by 7.2; and the truncation 
correction, the rounding correction, and the nonanalyticity correction vec- 
tors for the corresponding errors are denoted by ya'z, ya"’z, and ya'"z, 
respectively. 

10. Differential Equation for the Error.—Since z + cz is the desired 
solution of the system of equations, the following relationship holds, i.e., 


(12) (Zu + cu)’ = faults + C21, --+, S14 + Cora, f). 


For the intervals between points of computation, z also satisfies the 
equations (11), and, consequently, 


(13) (c2y)’ = fulti + car, -:+, S14 + Cora, t) — fale, --+, 214, 2). 


As 2 is now considered to be known, equation (13) is a differential equation 
on the vector, cz, which holds between points of computation. At points of 
computation, cz has a jump, yz, corresponding to the error made at this 
point. This jump condition and equation (13) determine precisely the cor- 
rection, cz. (See Illustration 2.) 

For practical computational reasons, equations (13) are made linear, i.e., 


r= > oh 

(14) (czy)’ = L 35, 2 

The coefficients 0f,/02, are regarded as constants over an interval which is 
large relative to the length of a step. Undoubtedly, this introduces an error 
into the error computation. However, consider the result of solving (14) as 
simply the zero order approximation to a solution of (13) by a Picarp itera- 
tion process. The use of (14) can be justified by estimating the corresponding 
solution of (13). By using higher terms in the expansion of the right-hand 
side of (13), one can readily justify the use of (14) for the purpose of attaining 
the real objective, the estimation of the size of cz,. 

One effect of the linearization of the differential equations is that the 
correction vector, cz, for the truncation error, rounding error, and non- 
analyticity error can be separately computed, provided one error is assumed 
to have a negligible effect on the others. This is clear for an estimate of the 
truncation error and nonanalyticity error. This is also justified for the 
rounding errors since only maxima and variances for these errors are 
computed. 

Thus, a truncation correction, c’z, can be defined as satisfying (14) be- 
tween computation points and having a jump, ya’z, at the a computation 
point. A similar definition holds for c’’z and c's, the rounding and non- 
analyticity corrections. The sum c’z + cs + c’"z satisfies (14) between 
computation points and has the required jump, ya%, at the a computation 
point. 

A function, c,'z, is introduced which is zero before the a computation 
point, has a jump, 7.’s, at this point, and satisfies (14) after this point (see 
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Illustration 3). The definition of c,’"2 and c,'""z is similar, and one can readily 
show that 


(15) c's = Ya la’2, o's = Ya la’'2, C'S.@ 2 abe th 


11. Error Effect Procedure.—In view of formula (15), the immediate 
problem is to compute c,’z. As was previously mentioned, intervals will be 
considered over which the equations (14) can be regarded as linear differ- 
ential equations with constant coefficients. The following equation can then 
be written in vector form, 


(16) (cz)’ = Hez. 

This equation is solved by finding the characteristic roots, \1, ---, Ai, of H 
and the corresponding characteristic vectors. Although zero is a double 
characteristic root, there are 14 characteristic vectors, x, ---, x“. The 
general solution is in the form, 

(17) y=>D kee, 

Now ¢,’z is the solution which at ¢ = t, has the value y,’z. Let X™, ---, X@® 
be the set of dual vectors to the characteristic vectors x, ---, x“, ie., 
(18) x XM = 643. 

Then 

(19) bs = : i Val BX Derilttady 


. 


Of course, the \,’s, x“’s, X’s assume complex values, but c,’z will be real 
since y,’z is real. A similar formula holds for c,’’2 and ¢,'"’z. 

In computing the \,’s, a further assumption was made which permitted 
one to factor the determinant |H — \| into two quintics, a quadratic, and 
\*. The quintics were solved by finding a real root by NEwrton’s Method, 
and the remaining quartics were solved by a method of Lin.? 

The process by which the characteristic vectors and their duals were 
obtained was somewhat complicated, since it was desirable to avoid four- 
teenth-order matrix operations. However, this was specific to the problem 
and will not be discussed here. 

12. Truncation Error.—In view of the above mechanism for evaluating 
the effect, c,’z, of an error, y.’2, one has the problem of obtaining ,’z in such 
a form that the summation process of (15) can be effectively carried out. 
For 7.'2 the difference is found between the open and closed integration 
values at every tenth step. This was expressed in analytic form. However, 
there is a more effective way of obtaining y,’. 

The f; for each 7 is given. Furthermore, by inspection it is seen that each 
fi can be expressed approximately in the form 


(20) fi = fP(t) + are cos (Bit + 1) + a2e~%2' cos (Bat + v2) + ase’, 


where f(é) is a slowly changing function, while a; + 18:, a2 + 782, and Xo 
are roots of the characteristic equations for H. The truncation correction 
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which is needed when an open integration is performed on the function, f, is 


tth 
(21) f fat — i h[SSf(t) — S9f(t — h) + 37f(¢ — 2h) — Of(t — 3h)]. 


This expression for the correction is linear, and if f is expressed as a linear 
combination of terms in the form Ae™* (where \ is complex), then for each 
such term the error is 





Y i < = rs = h 2h __ QpSdh 
(22) hAe~ [ vi 74 (55 — 59e” + 37e Se) |. 
= Ae™hT (Ah). 
Thus, in general, if 
N 
f=Lae™, 
i=1 
the truncation correction is 
N 
(23) yz = > hae T (Ah). 
i=1 


From the point of view of frequency analysis, T'o(Ak) is exceedingly 
interesting. For, if I'o(u) is expressed as a power series in u, we obtain 


(24) To(u) = .349u* + .716u5 + .789u® + .613u? + .375u® + .1914u* 
+ .084u° + --- (see footnote 3). 


This reflects the fact that the method of integration used was set up to 
integrate cubic polynomials perfectly. From the frequency standpoint, this 
means that there is a fourth-order contact at A = 0. On the other hand, this 
suggests an investigation of other integration formulae which would be 
precise at certain other frequencies particularly chosen for the problem at 
hand. Such an investigation has been initiated with some interesting results. 
For instance, an arbitrarily good approximation cannot be found over a 
connected Ah interval of length exceeding 2x no matter how many constants 
are available in the integration formula; but for a smaller interval such an 
approximation can be found. 

The above formula indicates that the contribution of the low frequency 
terms to 7’z; is negligible, and indeed only one term makes an effective 
contribution to ’z;. Consequently, one can obtain ¢,’z by (23) and (19). 
Furthermore, by a summation, based on a geometrical progression, it is 
possible to find c’z; using (15). 

Since the same frequency appears in both (23) and (19), it is noted that 
when the summation (15) is carried out, a term in the expansion of cz, is 
obtained in the form, 


(kit + ke)e~=¥ cos (Bit + ), 


while the other terms are exponentials without / in the coefficients. As far 
as this term is concerned, the total error builds up until the rather small 
exponential, a:, eventually cuts it down. This was the important term in 
the error analysis. 
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Thus even though the errors, 7.’z, oscillate in sign, their total effect 
builds up. The result is analogous to the solution of a differential equation, 


L(y) = g(x), 


where L(y) is a linear differential operator with constant coefficients, whose 
indicial equation has roots with negative real parts and in which g(x) con- 
tains a term corresponding to one of these roots. This phenomenon produces 
a resonance effect. 

This resonance phenomenon will occur whenever the original system of 
equations (1) and the variational equation (16) are close enough so that 
terms with frequencies associated with the variational equations appear in 
the f’s of (1) regarded as functions of the time. If practical, the method of 
integration should discriminate against these frequencies. 

As was mentioned previously, a stretch from ¢ = 2.10 to t = 2.90 was 
available over which there is a solution corresponding to an interval of 
h = .04 and one corresponding to h = .02. It is readily seen from the above 
that the error for the smaller interval should be approximately #; that for 
the larger interval. Now the theoretically computed estimate of the error can 
be compared with the difference between these solutions. The agreement 
between the two curves was very good in both size and shape. 

Since the error is a solution of a system of linear equations, the improve- 
ment in y.’2 obtained by halving the interval will be reflected in c,z. Thus, 
halving the interval is an excellent method of estimating the error when h is 
small enough to justify the linearization upon which (14) is based. Unfortu- 
nately, the cost in machine time is tripled by this procedure. 

13. Rounding Error.—To find the correction vector, +’’z, for the round- 
ing error, it is necessary to compare the correct calculation with the actual 
calculation. Naturally, this involves a study of the computation process used 
for each f. The correction, ya’’2;, for the 7 component is made up of a sum 
of a number of terms, a few of which are periodic, but most of which can be 
regarded as chance variables. Since the total effect of rounding error in this 
problem is not great compared with the other errors, it was felt that certain 
simplifications were justified. Hence, the y.’’2 was considered to be in 
the form 


(25) Ya''% = (ay€1", A2e2%, +++, Dis€r4"), 


where the a; is a maximum value for the component involved and the ¢,*’s 
are independent evenly-distributed chance variables taking on values be- 
tween —1 and 1. However, in the analysis, a;e/* could be replaced by a more 
precise expression for this component if it were necessary. 

It follows from (25), (9), and (15), that the total effect of the rounding 
error is in the form 


(26) c's = } oa" Ea, ;(t)e;*. 


The maximum can be estimated by proper summation methods; and it is 
also possible using SCHWARZz’s inequality and these same summation 
methods to obtain a bound for the variance of c’’z. 
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14. Points of Nonanalyticity.—In those cases in which an f; was non- 
analytic, it could be written in the form 


(27) f; = F; + H; sign (2; — A), 


where F; and H; are analytic and A is a constant. In general, i ~ j, and it 
was a relatively simple matter to compute both the correct integral for the 
term, H; sign (z; — A) and the value obtained from the method of inte- 
gration used. There is a discrepancy due to the fact that this method of 
integration assumes that the integrand is a polynomial while the term 
H; sign (2; — A) is in general discontinuous. However, this difficulty per- 
sists over four points, and the successive errors obtained in this way tend to 
cancel. Therefore, except for a temporary disturbance, the effect of an indi- 
vidual nonanalyticity is slight. Since there are relatively few of these, the 
total effect of the nonanalyticities is readily calculated by means of (19) 
and (15). 

15. General Planning Considerations.—Experience suggests that the 
preliminary analysis of a problem of this sort be based on the matrix, 
H = ((0f;/82.)), which, of course, can be evaluated initially. The character- 
istic roots of H indicate the frequencies present. The total error permitted 
can be divided into three not necessarily equal parts, one for truncation, 
one for rounding, and one for miscellaneous. In view of the resonance error, 
the total truncation error should be divided by the expected length of the 
computation to obtain the error permitted per unit interval, €, provided 
the situation is stable. If the open integration formula is used and if one 
, 4’, predominates among the characteristic roots of H, then h, the interval 
width to be used initially, can be obtained from 


(28) & = To(d’h). 


If more than one \ is important in H from this point of view, €o will have to 
be further subdivided and a portion allocated to each effective root. If some 
d has positive real parts, i.e., if the situation is unstable, it is usually per- 
missible to make a change in scale and a similar discussion applies. 

After having left the region around the initial value, # can be determined 
by comparing the result of an open and closed integration at specified inter- 
vals. In this connection, it should be pointed out that the quantity, T'.(Ah), 
which determines the error for a closed integration in the same way as Ip 
did for an open integration, is given by 


(29) — (.0264u* + .0341u5 + .0246u* + .0129u7 + .0055u! 
+ .0020u° + .0006u" + ---). 


The change in sign again indicates that the results of an open and closed 
integration span the true value, and it is apparent that the difference is an 
excellent estimate of the truncation error. 

16. Conclusion.—It seems clear from the above that a plan of solution 
should be based on considerations of over-all error rather than upon con- 
sideration of error per step. One such plan is given. This plan indicates a 
method for choosing h initially and also during the course of the solution. 
The ideas used can be supplemented by algebraic methods to estimate the 
total error when the solution process has been completed. 
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Our experience involved a stable system of differential equations, but the 
use of frequency analysis is justified in most cases, including unstable cases 
in which a scale change occurs. 

F. J. M. 

1 UK University, Computation Laboratory, Annals, v. 16, p. 176-187. [MTAC, 
v. 3, p. 437. 

2 SatH-NGE Ltn, “Numerical solution of complex roots of quartic equations,” Jn. Math. 
Phys., v. 26, 1947, p. 279-283. 

3 This formula assumes that the computed and hence available f is the function to be 
integrated. If one wishes to make an allowance for the distinction between this f and the 
correct f associated with the true solution, a difference equation must be solved. In case 
is a characteristic root of H, the corrected value Ip is obtained by dividing the above formula 

) by 
1 + .5000u + .1667u? + .0417u? — .0972u4 — .1424u5---. 


This correction is significant. 


RECENT MATHEMATICAL TABLES 


761[A].—R. Cousrat, “Calcul de V2, et réflexion sur une espérance mathé- 
matique,”’ Acad. Sci. Paris, Comptes Rendus, v. 230, 1950, p. 431-432. 


The first four terms of the binomial expansion of 


V2 = a(1 — 2x)-4 
where 
a=2—4&* 


and a@ is an approximation to V2, good to 333D, were used to obtain ¥2 to 
1032D. Besides this value the author gives the distribution of digits in the 


1033S values of V2 and 1/¥2. In the first 1000D in the V2, the digits 0-9 
have the following frequencies 


108, 98, 109, 82, 100, 104, 90, 104, 113, 92. 


Such a distribution has a chi-square of 8.38. The probability of such a value 
from a normal distribution is almost exactly 1/2. For 1/ V2 the probability 
is merely .05. [Compare MTAC, v. 4, p. 109-111]. 

The author ‘‘reflects’’ on the paradox that if one takes the product of 
the first 1033 digits of the decimal expansion of a real number x in the 
interval 0 < x < 1, the expected value of the product is (9/2) > 10%, 
whereas the probability that it is exactly zero is 1 — (9/10) > 1 — 10-*. 

D. H. L. 


762(C].—H. S. Unter, “A mathematician’s tribute to the state of Israel,” 
Scripta Mathematica, v. 14, 1949, p. 281-283. 
The author gives In 173 and In 5709 to 290D. 

763[D, H, L].—C. N. Daviss, ‘“‘The sedimentation and diffusion of small 
particles,’’ R. Soc. London, Proc., v. 200, 1949, p. 100-113. 
This paper contains a table of the first 16 positive roots of the equation 


2ax + tanx = 0 
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for a = 2, 1/2, 1/3, 1/6, 1/30. Values are to 5S, 6S and 7S. For a = 1/2 the 
equation is tan x = — x and the first 11 roots given in this case agree with 
those of PooLer, MTAC, v. 3, p. 496. 


764[D].—Tosnizo Matsumoto, “On Hayami’s turbulent tensor,” Kyoto, 
Imperial University, College of Science, Memoirs, v. 24A, no. 2, 1944, 
p. 63-72. 


On p. 66-69, is a table, computed by Hrros1 NAKAHETA, of the functions: 
z' cos $22, to 4D, 5; z! to 8D; cos 4x2 to 7D, for z = 0(.01)2. 
R. C. A. 


765[E].—H. G. Hopkins, “Elastic deformations of infinite strips,”” Cam- 
bridge Phil. Soc., Proc., v. 46, 1950, p. 164-181. 


The appendix (p. 181) is a 6D table of 


g1 = (16 cosh? x)/(55 cosh? x + 25x? + 9) 
and 
gs = (16 cosh’ x)/(39 cosh? x + 9x? + 25) 


for x = 0(.1)6. At x = 6 the functions are already near their limiting values 
16/55 and 16/39. 


766[F].—D. JARDEN & A. Katz, ‘“‘Additional page (477) to D. N. Lehmer’s 
Factor Table,” Riveon Lematematika, v. 3, 1949, p. 49 [English sum- 
mary p. 52]. 

This is the same as UMT 85[F], MTAC, v. 4, p. 29. 


767[F].—A. Katz, “Some more new factors of Fibonacci-numbers,”” Riveon 
Lematematika, v. 3, 1949, p. 14 [English summary, p. 54]. 


The author continues the factorization of the terms of the Fibonacci 
sequences U, and V, [see MTAC, v. 3, p. 299] giving the complete factoriza- 
tion of Uj,7 = 2-233-29717 - 135721 -673024656781, Vrs3= 151549 - 11899937029 
and Vios = 2-7-23-6263-103681-177962167367. The factors 128621 and 
119809 are given for Viog and Vi2s respectively. No further factors of the two 
series exist below 2-10° up to m = 128. 

ie es 


768[F].—G. PatamA, ‘‘Tabella delle posizioni iniziali relative al ‘Neocribrum’ 
di L. Poletti,’”’ Parma, Univ., Rivista Mat., v. 1, 1950, p. 85-98. 


The “Neocribrum” is a form of factor table devised by L. PoLett1 
[MTAC, v. 3, p. 532] and has for column headings the set of 6 numbers 


(S) 1, 7, 11, 13, 17, 19, 23, 29 


which are those prime to 30. Two cells of the table in the same column and 
adjacent lines correspond to numbers differing by 30. Once a prime p > 5 
appears as a factor of a number in any one column it reappears # lines farther 
down and continues to appear periodically. Hence to construct the factor 
table it is necessary to know the line number at which a given prime factor 
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will first appear. The purpose of the paper is to supply this information for 
all primes p for which 17 < » < 3547 (the primes 7, 11 and 13 are already 
printed into the otherwise blank forms of the sieve). Hence the table gives 
for each of the 491 primes the solutions x; of the congruences 


30(x; — 1) +m;=O0(modp) (¢ = 1(1)6) 


where m; are the numbers of (S). The awkward x — 1 is due to the fact that 
Poletti numbers his lines beginning with 1 rather than 0. The table gives 
also the least positive h for which h + 1001 is divisible by p. This is to enable 
the application of the results to successive ‘‘cycles”’ of the table. 

It is a little hard to see the need for publishing such a table. If one is 
going to construct a factor table of just this kind it is clearly indispensable. 
However, the table has other uses. The column headed 29 gives directly the 
value of 1/30 modulo # and the first column is a convenient list of primes. 

D. H. L. 


769[G].—F. N. Davip & M. G. KENDALL, ‘“‘Tables of symmetric functions— 
part I,” Biometrika, v. 36, 1949, p. 431-449. 


The authors present tables of coefficients in the linear representation of 
the monomial symmetric function in terms of products of sums of like 
powers, and conversely the coefficients in the expansion of the latter in terms 
of the former. The tables extend to symmetric functions of weights <12. 
Thus from the two tables of weight 4 we read, for example, that 


22 aBy = 22 ot — 2(D o)(L a) — (2 a”)? + (L @*)(Z a)? 
and that 


(X a*%)(L a)? = Dot +22 a8 + 22 af? + 2D By. 


In order to avoid fractions the monomial symmetric functions are multi- 
plied by the product of factorials of the exponents in the partition of the 
weight represented by the monomiai. The resulting functions are called 
- “augmented.” Thus >> a°6*y*e is multiplied by 3! 2! 

The authors appear to be unaware of the tables of SUKHATME, ZIAUD- 
Din, KErawaLa and Kerawata & HanaFi [MTAC, v. 3, p. 24]. The first 
of these gives the same information for weights <8. The others give only 
half as much information (i.e., the coefficients in the expression of the 
monomials in terms of sums of like powers) for weights 9, 10 and 11. The 
table of weight 12 under review is completely new. 

The various checking procedures used by the authors should be sufficient 
to produce a set of tables completely free from error. Nevertheless it would 
have been wise to collate the tables with those previously given which are 
known to contain errata [MTAC, v. 3, p. 24]. The application of the tables 
to problems in statistics is treated briefly. 

The tables are arranged beautifully in lexicographical order, the two 
triangular halves fitting into a perfect rectangular layout. The type, though 
rather small, is quite clear. 

D. H. L. 
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770(G].—J. A. Topp, ‘The characters of a collineation group in five dimen- 
sions,” R. Soc. London, Proc., v. 200, 1949, p. 320-336, insert between 
p. 336-337. 


The paper gives a table of characters of a primitive group of collineations 
of order 6531840 = 2*-3*-5-7 in five space and a table of characters of a 
subgroup of index 2. The group is described in another paper.’ 


1J. A. Topp, “The invariants of a finite collineation group in five dimensions,” Cam- 
bridge Phil. Soc., Proc., v. 46, 1950, p. 73-90. 


771{I]|.—AnveErs Rez, “On quadrature formulae,’”’ Cambridge Phil. Soc., 
Proc., v. 46, 1950, p. 119-126. 


The author considers the usual quadrature formula for an integral with 
weight function w(x): 


b n 
f w(a)fiz)de = ¥ pefles) + Re 


As is well known, the best choice of x; is the set given by the Gauss’ method 
corresponding to w(x). However, these values are not rational and the com- 
puter is faced with the task of interpolating to find f(x,). If we choose for x; 
the Gaussian values rounded off to 2D and modify the p; accordingly we 
obtain a formula which is theoretically easier to use and which retains 
nearly the full force of the Gauss approximation. 

Six tables of the x;, ; are given for the six weight functions 


, witi-s, #(i-s) (1 — x)X1 + 2)-, 
e and ax exp (—x’). 


The order m extends as far as 10, 8, 8, 5, 5, and 8 in these respective cases. 

For the first four weight functions the interval (a, 5) of integration is 
either (0, 1) or (—1, 1) while in the last two cases it is (0, ©) and (— @, «) 
respectively. 

The x; are given to 2D, as mentioned above, and the coefficients p; to 7D. 

Five examples are given comparing the present method with the strict 
Gauss formula and WeEDDLE’s rule. As might be expected, the modified 
Gauss method is only slightly inferior and considerably easier to use. It is 
much superior to Weddle’s rule. 

In “actual practice’ the above avoidance of interpolation difficulties may 
be, in some cases, only apparent. In fact if the function to be integrated is 
observed or tabulated on an interval (a, 5) different from (0, 1) or (—1, 1), 
such as (0, 27), then a transformation of variable is required which may 
introduce interpolation after all. 

D. H. L. 


772(T|.— HERBERT E. Sauzer, “‘Formulas for complex cartesian interpolation 
of higher degree,” Jn. Math. Phys., v. 28, 1949, p. 200-203. 


The purpose of the present table is to provide the coefficients for a 
Lagrangean interpolation polynomial adaptable to interpolation over a 
square grid in the complex plane. The polynomial assumes the form 


) a,F (zx) 
f (P ) - Lia ’ 
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where 3% = k are the points of grid-configuration, a, = A;/(P — k), and 
P = p + qi is the variable point. For example, in the case of three-point 
interpolation the configuration consists of the three points 0, 1, i and the 
values of A; are correspondingly Ay = —i, A; = $(1 + 1), A; = $(—1 +2). 
From this we find that >> a reduces to 1/[P(P — 1)(P — i)] and the 
interpolation polynomial becomes 


S(P) = —(P — 1)(P — Ofot+ 41 + OPP — Of + 2(-14+0P(P — 1)f.. 


Beginning with five-point interpolation, the author gives alternative 
configurations ‘‘which have the property of making the location of P (or z) 
considerably more central with regard to the points k (or z), and hence 
would be expected to yield greater accuracy. The number of these more 
central configurations which are given are: for five-point—one, for six- and 
seven-point—two, for eight- and nine-point—three. Thus, the user has a 
great latitude of choice in available formulas for complex interpolation, for 
checking that interpolation, and for central choice of the argument P.” 

This paper is a continuation of the original work of A. N. Lowan and 
the author! [MTAC, v. 1, p. 358-359] and a paper by the author.? The work 
makes use of a simplification introduced by W. J. TayLor® for real Lagran- 
gean interpolation. 

H. T. Davis 
Northwestern University 
Evanston, I!linois 


1A. N. Lowan & H. E. Sauzer, ‘Coefficients for interpolation within a square grid in 
the complex lane,” Jn. Math. Phys., v. 23, 1944, p. 156-166. 

2H. ALZER, “Coefficients for complex quartic, quintic and sextic interpolation 
within a square grid,” Jn. Math. Phys., v. 27, 1948, p. 136-156. 


. J. Taytor, “Method of Lagrangian "curvilinear interpolation,” NBS, Jn. Res., v. 
35, 1945, p. 151- 155. 


773(K].—D. J. Gres & J. N. Berretront, “AOQL single sampling plans 
from a single chart and table,”” Amer. Stat. Assn., Jn. v. 44, 1949, 
p. 62-76. 


As the title states, a single chart and table are given in the paper to find 
the AOQL (Average Outgoing Quality Limit) single sampling inspection 
plan which will yield a “practical” minimum amount of total inspection. 
Total inspection is used here to mean the combined amount resulting from the 
single sample inspected from each and every lot and the 100% screening in- 
spection of lots rejected under the single sampling plan. Given an AOQL value 
in % it is desired to maintain and the lot size, N, Chart I of the paper (p. 66) 
is entered to find the acceptance number, c. Table II of the paper (p. 65) 
is then entered along with the AOQL value to find the sample size m for the 
single sampling plan (c, m) with c = 0(1)12, AOQL = .1,.25(.25)1(.5)5(1)10. 

The present paper and associated tables appear to eliminate the necessity 
of knowing the process average accurately, at least for many practical 
situations. 

F. E. Grusss 
Ballistic Research Laboratories 
Aberdeen Proving Ground, Md, 
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774[K].—D. J. Fuyney, “The truncated binomial distribution,” Annals of 
Eugenics, v. 14, 1949, p. 319-328. 


A number s of observations of binomial random variables, all with the 
same probability p of success, are obtained under circumstances which make 
it possible to observe only values greater than zero. The problem is to obtain 
the maximum likelihood estimate of ». Referring to his earlier paper’ for 
derivations, the author presents an iterative technique for solving the likeli- 
hood equation which is made simple by giving a 3D table of the weights of 
single observations and of the bias in weighted scores for s = 2(1)20 and 
p = .01(.01).05(.05).95. A similar technique is described for the case of 
doubly truncated binomial distributions in which neither the number 
of ‘“‘no successes” nor the number of “no failures’ can be observed, and 
a similar table of weights and biases, both to 3D, for s = 3(1)20 and 
bp = .01(.01).05(.05).50 is given. 


University of Washington 
Seattle, Washington 


1D. J. Finney, “The estimation of the aed of recombinations. I. Matings of 
known phase,”’ Jn. Genetics, v. 49, 1949, p. 159-175. 


Z. W. BiIRNBAUM 


775[K].—_N. L. Jonnson, ‘Systems of frequency curves generated by 
methods of translation,” Biometrika, v. 36, 1949, p. 149-176. 


The author investigates the properties of probability functions, »,, where 
z = y + 6In f(y) is distributed normally with unit variance, f(y) = 0. If 
f(y) =y we obtain the well-known log-normal system, S_; if f(y) = y/(1 — 9), 
0 < y <1, a new system Sz; and if f(y) = y+ Vy? +1, — @ <y< @, 
the new system Sy. On page 157 is a chart in terms of the PEARSON measures 
of skewness and kurtosis, 8; and 82, showing the regions of S,, Sg, and Sy. 
On page 164 is a nomogram which gives 6 and 7/4 in terms of 8; and Bs, 
0 < 6: < 1.3, 3 < B. < 5 for the system Sy. Table 8, page 174, gives the 
values of yi’, o, 81, Be for 6 = .5, 1, 2, and y = 0(.5)2.5 for Sg. The author 
applies his results to the graduation of observed frequency distributions and 
to the normalization of skewed distributions. “ee 


Hughes Aircraft Company 
Culver City, California 


776[K].—K. R. Narr, “A further note on the mean deviation from the 
median,”’ Biometrika, v. 36, 1949, p. 234-235. 


Values of the coefficient of variation (2S) are given for the mean devia- 
tion from the mean and the mean deviation from the median for samples of 
size 2(1)10 from a normal population. They are the same within 1 figure in 
the second place. This implies that the mean deviation from the median is 
just as precise an estimate of dispersion as the mean deviation from the 
mean for samples from a normal distribution. 

Actually there are more precise linear estimates of dispersion than either 
of these. Some examples are given in a text by Drxon and the reviewer.! 


University of Oregon F. J. MAssEy 
Eugene, Oregon 


1W. J. Drxon & F. J. Massey, Introduction to Statistical Analysis, lithograph ed. 
Eugene, 1949. 
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777[K].—E. G. Oups, ‘‘The 5% significance levels for sums of squares of 
rank differences and a correction,’ Annals Math. Siat., v. 20, 1949, 
p. 117-118. 


In 1938 the author published! a series of tables having to do with the 
distribution of the rank correlation coefficient. Table V gave pairs of 
values between which }> d? (d; being the rank difference for the 7“ indi- 
vidual) has a probability, P, of being included under the hypothesis that 
x @d = (n* — n)/6. This hypothesis is equivalent to the null hypothesis that 
the rank correlation coefficient is zero. The table gave pairs of values for 
n = 11 tom = 30 inclusive and for P = .99, .98, .96, .90 and .80. The table 
in the present paper extends the original table by giving pairs of 1D values 
corresponding to P = .95. 

The correction is for the formula printed in the earlier publication for 
the variance of the normal deviate, x, used in calculating the values of the 
table. 

CarL F. Kossack 
Purdue University 
Lafayette, Indiana 


1E. G. Oxps, “Distribution of the sums of squares of rank differences for small num- 
bers of individuals,’’ Annals Math. Stat., v. 9, 1938, p. 133-148. 


778[K].—FRANCES SWINEFORD, ‘Further notes on differences between per- 
centages,’ Psychometrika, v. 14, 1949, p. 183-187. 


The two tables are designed to determine the least common size, N, for 
each of two samples for testing the hypothesis that the difference in the 
two population proportions, p: — p2, is at least d;. It is assumed that the 
sample proportions, ~,’ and ,’, are distributed normally and, therefore, that 
the appropriate test is a one-tailed test of the hypothesis that p; — p2 = d,. 
Then N = 5.4119(p191 + poq2)/(do — ds)? at the 1% point and approxi- 
mately half as much (1/2.0003) at the 5% point, where dy = p,' — p,’. The 
tables give N for the 1% points only. 

Table 1 gives N’ = 10.8238pq/(dy — d,)*, where p = $(p1’ + p2’), to 
OD for p = .10(.05).90 and |d» — d;| = .050(.002).080(.005).135. Table 2 
gives the correction terms (191 + $292)/2pq to 3D(3S) for p = .10(.05).90 
and d; = .10(.05).50. 

Leo Katz 
Michigan State College 
Lansing, Michigan 


779[K].—J. E. Watsu, ‘‘On the power function of the ‘best’ ¢-test solution 
of the Behrens-Fisher problem,” Annals Math. Stat., v. 20, 1949, p. 
616-618. 


Let m sample values be drawn from N(a,, ¢;*) and m sample values from 
N(a2, 027), m < n, (where N(a, o”) represents a normal probability function 
with mean a and variance o*), then SCHEFFE has shown that a ‘‘?’”’ test solu- 
tion of the BEHRENS-FISHER problem, ¢;?/¢;? not known, using in the numer- 
ator the difference of sample means, and the denominator based on the 
square root of a function of the sample values which has a x*-distribution 
with m — 1 degrees of freedom has certain optimum properties. The purpose 
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of the note is to compare the power function of this ¢ test with the power 
function of the correspondingly most powerful test for the case in which the 
ratio ¢;7/¢;? is known, for one-sided and two-sided symmetrical tests. This is 
done by finding the power efficiency. Two 3D tables are given for the power 
efficiency of Scheffé’s test, o:7/o7 not known, against the test ¢;*/¢,7 known, 


on page 617 where a = significance level, for a = .05, m = 4, 6, 10, 15, 20, 
30, 50, 100, « and the same range for m; also a = .01, m = 6, 8, 10, 15, 20, 
30, 50, 100, ©, same range for n. 

L. A. AROIAN 


Hughes Aircraft Company 
Culver City, California 


780[K].—J. E. Watsu, “On the ‘information’ lost by using a t-test when the 
population variance is known,” Amer. Stat. Assn., Jn., v. 44, 1949, 
p. 122-125. 


Table 1 gives the approximate number of sample values ‘“wasted”’ if, 
when the population variance is known, one uses a t-test (estimating variance 
from sample) in place of the appropriate normal deviate test, when one is 
testing whether the population mean differs from a given cofstant value. 
5%, 2.5%, 1% and .5% significance levels are tabulated to 2S for both the 
one-sided and symmetrical test. Sample values wasted is defined in terms 


of equal power functions. a 


Purdue University 
Lafayette, Indiana 


781[K].—J. E. Watsu, “Applications of some significance tests for the 
median which are valid under very general conditions,” Amer. Stat. 
Assn., Jn., v. 44, 1949, p. 342-355. 


The tests considered are valid for samples of m from one or more uni- 
verses with a common median, which are symmetric and have continuous 
cumulative distributions. Table I lists one to five such tests for each 
nm = 4(1)15 with their approximate one-sided and symmetric significance 
levels to 3D and their efficiencies to nearest .5 on the assumption that the 
universes sampled are normal. Table II lists further tests for » = 4(1)9 for 
which the bounds of these significance levels are given, as well as their sig- 
nificance levels and efficiencies to same precision as Table I for samples 
from normal. Table I was also published elsewhere.' 

University of Mississippi T. A. Bocxunstare 
University, Miss. 

1J. E. Watsu, “Some significance tests for the median which are valid under very 

general conditions,”” Annals Math. Stat., v. 20, 1949, p. 64-81. 


782[L].—S. CHANDRASEKHAR, “‘On Heisenberg’s elementary theory of turbu- 
lence,” R. Soc. London, Proc., v. 200A, 1949, p. 20-23. 


The function f(x) satisfies a certain nonlinear differential equation. 
Introducing the new variables 


y= [sora and ¢ = 2f(2), 
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and denoting differentiation with respect to y by primes, the differential 
equation becomes 


gig’ + 2y(4 + g’) + 2g8(4 — g’) — 8g = 0. 


gly) = 4y + yi(a + (4/3) ny) + -:- 


where a is an arbitrary constant. Starting with this value, the differential 
equation was integrated numerically for a = 1.8104739, 1.81, 1.75, 1.5, 1.0, 
.5, 0. Four decimal tables of f(x) are given for various ranges of x and varying 
intervals. 


For small y, 


A. E. 


783(L).—L. HowartTH, ‘‘Rayleigh’s problem for a semi-infinite plate,’’ Cam- 
bridge Phil. Soc., Proc., v. 46, 1950, p. 127-140. 


The work leads to the expression 


ome uW 1 why x (1/8) v? dv 
oa (vt)! aw + aa I, “ K,(4v’) Pe 


and the notation rg = (xvt)~*yW is used. Table 1 gives values to 3 decimal 
places of (W?t/v)!r/(pW?) and of (uW) Jy’ (7 — re)dr for R = (vt) = .01, 
05, .1, .2(.2)2(.5)4 and o. 

A. E. 


784{[L].—_M. K. Krocpau., “On broadening of hydrogen lines in stellar 
spectra—I,” Astrophys. Jn., v. 110, 1949, p. 355-374. 


The function 
a(§ 2 (3) 2 aa(3) 2 12 2 
~qs\2) St s\2) Ms 30\2) 262s 
42° dy sf si dy 2(3)'[- dy 
et. Tee Mile Rie MND de: oy Vy 
is tabulated to 3D for z = 0(.1)1(.2)3(.5)5, and to 2D for z = 6(1)10. The 
computations were apparently based on tables contained in JAHNKE-EMDE. 


[1933 ed. p. 6-9, 34-35]. 
A.E. 


785[L].—MERsT, ‘Untersuchung zur Arbeit von H. G. Kiissner ‘Lésungen 
der klassischen Wellengleichung fiir bewegte Quellen’ ’’ Aerodynamische 
Versuchsanstalt, Géttingen e. V., Abteilung J 06, Bericht B 44/J/31, 
ZWB/AVA/Re/44/J/31, ZWB 10514, 1944. English translation ““Wave 
propagation from moving sources” by O. W. Leibiger Research Labora- 
tories, Petersburg, N. Y., ATI No. 32413, 8-8-701, 1948. 14 p. 


The function 


1 2 
Iimn(a, o) = xf (1 — 20 cos x)~* exp {imn[— x — a(1 — 20 cos x)*]}dx 














tial 


tial 
1.0, 
ring 


am- 


imal 


01, 


ellar 


ingen 
ische 
j/31, 


Nave 
bora- 


idx 





RECENT MATHEMATICAL TABLES 153 


occurs in an investigation! of sound propagation from a source moving in a 
circle. In the present paper the following notations are used: 


B = imna, A = (1 — 2¢ cos x)}, y = Peooex, 
F(y) = A-e*4, —s f(a) = #4 


First the first 8 derivatives of F(y) are given in terms of 8 and A, then the 
values, in terms of 8, of these derivatives when y = 1. From these the first 
7 coefficients in the asymptotic expansion 


limn(a, ¢) = 4 > b,Jmn(vmnac) 

v=l 
are computed. Next, the first 6 derivatives of f(c) are given in terms of 8 
and A, and the values of these derivatives for ¢ = 0 in terms of 8 and x. 
The Maclaurin series of f(c) leads to an expansion of h in powers of ¢ which 
is said to show satisfactory convergence, at any rate for small m, for all 
values of a and ¢ which are of importance; in fact it is stated that a few 
terms of the series suffice for numerical computation. Certain polynomials 
which occur in the first 7 coefficients of the expansion of h,, in powers of ¢ 
are given explicitly. There are 14 such polynomials, 7 to be used for even m 
and 7 for odd m. Table I of the appendix gives the coefficients of these 
polynomials, and Table II certain auxiliary quantities which can be used for 
a rapid computation of the coefficients. 

A. E. 


1H. G. Ktssner, “Lisungen der klassischen Wellengleichung fiir bewegte Quellen,” 
Z. angew. Math. Mech., v. 24, 1944, p. 243-250 


786[L].—F. W. J. OLveR, “Transformation of certain series occurring in 
aerodynamic interference calculations,’ Quart. Jn. Mech. Appl. Math., 
v. 2, 1949, p. 452-457. 


by = Sauts © (—1)"Ke(2musn) + 4X (—1) aS 


n=l n 
is given in Table 1 (p. 456) to 6 decimal places for wu = .5,s = 1(1)6; 4 = 1, 
s = 1, 2, 3; and » = 2, s = 1, 2. 


Table 2 (p. 457) is an auxiliary table giving values to various degrees of 
accuracy of K;(mx) and IJ,;(mx) for 7 = 0,1 and m = 1(1)7. 


787(L].—F. Rrecets, “Formeln und Tabellen fiir ein in der raumlichen 
Potentialtheorie auftretendes elliptisches Integral,” Archiv d. Math., v. 
2, 1949-50, p. 117-125. 


The integral in question is 
i ~—cos 2n8 
GH) = (—1)) f G—Fsinto™ 


and can be expressed in terms of complete elliptic integrals of the first and 
second kinds in the form 


(1 — B) kf E(k*) — gnK(k*) J, 
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where f, and g, are polynomials of degree m in k”. The exact values of the 
coefficients of these polynomials are given (as common fractions) for 
n = 0(1)7. 

(2/x)G,(k*) can be expanded as a power series in k?, and 5 decimal values 
of the coefficients of this power series up to the coefficient of k® are given 
for m = 0(1)7. 

k”G,(k*) can also be expanded in the form: power series in k” + In (4/k’) 
times a power series in k”. The coefficients occurring in this expansion up 
to (and including) the terms in k’” are tabulated to 5 or 6 significant places 
for m = 0(1)7. 

The principal tables are 4 decimal tables of k” and of k”G,(k*) for 
n = 0(1)3, k? = .00(.01).99 and k? = .900(.001).999. 

The computations were carried out on the Hollerith equipment of the 
Max Planck Gesellschaft of Géttingen. The original computations used the 
interval .001 throughout and included also G, to G;. The unpublished parts 
of the tables are available through the author. 

A. E. 


788[L].—G. ScCHWEIKERT, “‘Zur Theorie des Gasdrucks gegen eine bewegte 
Wand,” Zeit. angew. Math. Mech., v. 29, 1949, p. 289-300. 


Let 


F(x) = 2a f exp (—/)dt 
and 
G(x) = wix exp (x?) F(x). 
The article contains tables to 5D of 
F(x), 2x? F(x), (2 — F)(1 + 2x?) + 2x-'x exp (—x’?) 

and 

[1 + 2x2(1 — 6) ]F(x) — 2x-ix exp (—x’) 
for 

6 = 0(.1).4, 1 and for x = 0(.1)1.5, 2, 2.5 


and tables to 8D of log G(x), 1/G(x) and 


1 + 3x? — 1/G(x) 
for x = .05, .1(.1)2(.5)4.5. 


789[L].—GEOFFREY TAYLOR, “The formation of a blast wave by a very 
intense explosion. I. Theoretical discussion.”” R. Soc. London, Proc., 
v. 201A, 1950, p. 159-174. “II. The atomic explosion of 1945,” ibid, 
p. 175-186. 


In the course of this work there appear the functions f, ¢, ¥ of » which 
satisfy the system of nonlinear differential equations 


VY _ o + 26/n 
v n—@¢ 


staf +h (-atef-of =0 


in which y is a constant. 
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Table i of Part I (p. 164) contains 3 decimal values of f, ¢, ~ for 
n = 1(—.02).5. These values were computed by step-by-step numerical 
integration of the differential equations with y = 1.4 and the initial values 
f(1) = 1.167, 6(1) = .833, ¥(1) = 6.000. 

Approximate formulas for the functions in question are developed on 
p. 165. 

Table 2 of Part I (p. 166) contains 3 decimal values of f, ¢ and 2 decimal 
values of y for 7 = 1, .95, .9, .8, .7, .5, 0. These values were computed from 
the approximate formulas with y = 1.666, f(1) = 1.250, $(1) = .750, 
¥(1) = 4.000. 

Table 2 of Part II (p. 178) contains 3 decimal values of f, ¢, y for 
n = 1(—.02).9(—.05).4 obtained by approximate calculation for y = 1.30, 
f(1) = 1.130, (1) = .869, ¥(1) = 7.667. Some values of the temperature 
are added. 

There are also other tables. 

A. E. 


790[L].—A. vAN WIJNGAARDEN & W. L. ScHEEN, “Table of Fresnel inte- 
grals,” Akademie van Wetenschappen, Amsterdam, Afd. Natuurkunde, 
Verhandelingen, eerste sectie, v. 19, no. 4, 1949, 26 p. (Report R49 of 
the Computation Department of the Mathematical Centre at Amster- 
dam.) Price 2.50 guilders. 


This is a table of 
C(u) = f cos $xf*dt and S(u) = f sin }xf*dt 
0 0 
for u = [0(.01)20; 5D] with modified second differences. 


(1) C(u) = x Carpet, Crys = (—1)*(9)*/[ (2k) (4k + 1)] 


(2) S(u) = Y Saysu***, = Sunes = (—1)*(54)***/[(2k + 1)1(4k + 3)]. 
k=O 

Tables are given for k = 0(1)22 or 21, of Cas:lCi = 1, Cs = — 2.467401- 

100272340-10-1, ---, Cog = 2-10-], and of Sais[. Ss = 5.235987755982989 


-107, ---, Sx = — 5 X 10-7]. 
Asymptotic expansions are 


(3) Clu) ~ $+ sin $x? > Varga OD — cos $a? FS ou 
k=0 ku 


(4) S(u) ~ $ — cos xu? > Veep ut) — sin $4? ye Can st) 
k=0 ko 

where 

(S) Yas = (—1)*2(4k) !(24)-O*+Y /(2k)! 

(6) Curss = (—1)*2(4k + 2)!(20)-@**”)/(2k + 1)! 
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There are tables of the coefficients (5) and (6) for k = 0(1)14 and 
k = 0(1)13 respectively. 

In computing the main table 9D values were calculated by means of (1) 
and (2) for « = 0(.5)2.5 and the same from (3), (4) for # = 2.5(.5)12. In 
order to get this accuracy by means of asymptotic series for values of u as 
low as 2.5, the technique described by GooDWIN & STATON!' was employed. 
The next step was to prepare preliminary 7D tables of C(u) and S(u) for 
u = 0(.01)12 by numerical integration of 5D values of cos }xf* and sin }rf* 
with an interval h = .01. After most elaborate treatments and testings, the 
authors were led to their rounded-off 5D values, which are ‘‘guaranteed.” 

The values of the table for u = 12(.01)20 were calculated by means of 
the asymptotic series and checked by complete duplication of the computa- 
tion. Modified second differences, *®* = 5? — 0.1845 were computed from 
the 7D tables, then rounded off to 5D and checked by differencing. 

Linear interpolation will yield no larger error than 4 X 10 X uw. Full 
profit of the accuracy of the table is obtained by the use of EVERETT’s 
formula up to second modified differences: 


f(xo + ph) = (1 — p)fo + Ofi + Ec%ic™* + E175;*. 


A small table of the interpolation polynomials E,? and E;* is given. 
Extracts from text. 


We have recently referred in MTAC to several other tables of C(u) and 
S(u): C. M. Sparrow (v. 3, p. 479), for uw = 0(.005)8; 4D; D. L. ARENBURG 
& D. Levin (v. 3, p. 479), for « = 0(.1)20, and wu = 8(.02)16; U. S. Navy, 
Res. LaB., Boston (v. 3, p. 417), for wu = 0(.1)20; 4D or 4S; R. T. Brrce 
(v. 4, p. 30), for uw = 0(.05)12.05; 4D. 

A comparison of the first 100 values of C(u) and S(u) of tables under 
review with the corresponding entries in the Sparrow tables indicated the 
following 29 apparent unit-errors in the fourth places of Sparrow: S(.15), 
C(.23), S(.31), C(.33), S(.33), C(.45), C(.47), C(.49), S(.57), S(.59), S(.63), 
S(.65), S(.69), C(.72), C(.73), C(.74), C(.75), C(.76), S(.77), C(.79), S(.79), 
C(.84), S(.85), C(.91), C(.93), S(.93), S(.94), C(.95), C(.98). There are also 
6 2-unit errors at C(.77), C(.78), C(.80), C(.81), C(.82), C(.83). These results 
suggest that tables based on Sparrow are likely to be not without error. 

Wijngaarden & Scheen make no reference to any earlier table of Fresnel 
integrals. A misprint for C(4.95), 0.45404 has been corrected by hand to 
0.54504. 

R. C. A, 


1E. T. Goopwin & J. Staton, “Table of {{° e“du/(u + )," Quart. Jn. Mech. Appl. 
Math., Oxford, v. 1, 1948, p. 319-326. See MTAC, v. 3, p. 483 


791[L].—D. V. ZaGREBIN, ‘‘K voprosu o tochnosti formuly Stoksa”’ [Con- 
cerning the accuracy of Stokes’ formula], Akad. Nauk. SSSR, Inst. 
Teoret. Astr., Bull., v. 4, no. 3 (56), 1949, p. 134-141. 


Table 1, p. 137, is a 3D table of six very special functions which are 
combinations of trigonometric and logarithmic functions with complete 
elliptic integrals. There are also tables to 3,4D of 15 definite integrals of 
products of these functions by even powers of the sine and cosine. 
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792[M].—Roy C. SPENCER, PAULINE AusTIN, ELIZABETH CHISHOLM, ELLEN 
Fine, & JEANE SCHWARTZ, Tables of Fourier Transforms of Fourier Series, 
Power Series, and Polynomials. Report S-58, July 10, 1945. Radiation 
Laboratory, Massachusetts Institute of Technology, Cambridge, Mass., 
ii, 29 p. 21.4 X 27.9 cm. 


The tables, occupying p. 6-27, are as follows: 


. I, p. 6-11, Fourier transform of a constant, go(¢?) = ¢~' sin ¢, g,’, its 
derivative Dgo, (Dgo)?, ¢ = [0(5°)1080°; 8D]. 

. II, p. 12-13, D*go for » = 1(1)8, @ = [0(30°)1080°; 8D). 

. III, p. 14, Fourier transforms T of cos mxx and sin mxx for n = 1(1)3. 
Tables of T cos nxx, iT sin nxx, n = 1(1)3, ¢ = [0(30°)720°; 7D]. 

. IV, p. 15, Fourier transforms of cos $nxx, sin 4nxx, n = 1, 3. Tables of 
T cos 4nxx, iT sin nxx, n = 1, 3, ¢ = [0(30°)720°; 7D]. 

. V, p. 16, T(1 — x*)*, m = 1(1)3, @ = [0(30°)720°; 7D]. 

. VI, p. 17, Fourier transforms, 7P,(x) = 1%j,(¢), m = 2(1)4. Tables of 
jn(o), m = 2(1)4, @ = [0(30°)720°; 7D]. 

. VII, p. 18-25, Fourier transforms of i*x", go, D*go, m = 1(1)8, 
@ = [0(.1)20; 8D]. 

. VIII, p. 26-27, Fourier transforms of (1 — x*)*, (1 + D*)"go, n 
@ = [0(.1)10; 8D]. 
There are a number of disagreements with values given in the tables of 

RMT 726 [MTAC, v. 4, p. 80-81]. 


3 8.459 8 FA A 


1(1)4, 


- ©. A. 


793[U].—PieERRE Hucon, Nouvelles Tables pour le calcul de la droite de hauteur 
a partir du point estime, Paris, Girard, Barrére and Thomais, 1947, xiv, 
92 p. 15.0 X 21.2 cm. + 1 chart 25 X 42.4 cm. 


The tables are two in number; the principal table was designed for use 
on shipboard in calculating the altitude for the dead reckoning position by 
logarithms and haversines. It consists of 90 pages of five-place values of the 
natural haversine, log haversine, and log cohaversine, with argument 
1°(1’)179°; each page contains an interpolation table for each of the three 
functions for tenths of a minute of arc. Preceding the principal table are two 
pages of corrections for refraction, height of eye and semidiameter to be 
applied to observed altitudes of the lower limb of the sun and of stars. The 
arguments are observed altitude 6°(1°)16°(2°)20°(5°)50°(10°)90° and height 
of eye 3(1)10(2)26 meters. 

The foreword and explanation are presented first in French and again in 
English. The formulas used are: 


hav (90° —h) = A+B 
log A = log cohav ¢t + log hav (d — L) 
log B = log hav ¢ + log cohav (d + L) 


where #, d, and h are the local hour angle, declination and altitude of the 
celestial body and L is the dead reckoning latitude. Since the haversine and 
the cohaversine are always positive and between zero and one, the rule of 
signs is no longer needed. Also since corresponding values of the natural 
haversine and log haversine are given side by side, no separate table of 
logarithms of numbers is given. 
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The azimuth is determined by the use of a nomogram folded inside the 
back cover; it is based on the formula: 


cosh cos Z = sin (d — L) cohav# + sin (d + L) havi 


where Z is the azimuth angle of the celestial body. It is intended that the 
azimuth shall be determined only to the nearest degree which is generally 
adequate for ordinary navigational purposes. 

The author suggests that the accuracy of his method as compared to the 
classical FriocourT method is as follows: 


h = 60° kh = 75° h = 84° 
Friocourt +0/5 +0/9 +2/2 
Hugon +0/4 +07 +1!7 


and hence the claim is for greater speed and ease of use rather than greater 
accuracy with the same number of decimals. 

The printing of the tables is rather poor on the whole, but a part of the 
trouble may be blamed on the quality of the paper which is mediocre. It is 
to be hoped that the proofreading of the tables has been done with greater 
care than that of the foreword and explanation. In the English explanation, 
X — Y = 1 should obviously be X + Y = 1, and in the expression for X 
in both the French and English explanations, cohaversine (D + ¢) should 
be cohaversine P. 

CHARLES H. SMILEY 
Brown University 
Providence, R. I. 


MATHEMATICAL TABLES—ERRATA 


In this issue references have been made to errata in RMT 790 (Wijn- 
gaarden & Scheen), 792 (Spencer et al.), 793 (Hugon), and Note 118. 


173.—Polnoe Sobranie Sochinentt P. L. Chebysheva [Complete Collection 
of Works by P. L. Chebyshev]. Volume 1, Teorita Chisel [Theory of 
Numbers], Moscow and Leningrad, Academy of Sciences, 1946. 342 
p. + portrait frontispiece. 15 X 23 cm. 20 roubles paper; 23 roubles 
bound. Edition (Second, stereotyped) of 3,000 copies. 


Previous editions of this volume have been reviewed in MTAC, v. 1, 
p. 440-441. The present volume not only reproduces the errata of the 1944 
edition but adds many new misprints both in the text and in the tables, 
p. 311-339. These latter are as follows: 


page 
311 line —20 for 2372 read 2237 
314,p=13, N=12 read I=6 
p= 19, insert N= 1 
p = 23, insert 4 between 17 and 5 
317, p = 61, WN table, for line 1 read 1 10 39 24 57 21 27 26 16 38 
p’=67, I=47 for N=38 read N=18 
318,p=71, N=16 for I=15 read I =22 
N=26 for I=22 read I =15 
319, p = 89, insert the primitive root 35 
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he page 

320, p = 103, [= 99 read N=3! 
322,p = 127, I= 116 for N=31 read N=T71 
323,p= 131, 1=35 for N=76 read N=79 
he 324,p=139,N=50 for I=26 read I =25 
325,p=151,N=12 for I=121 read I = 131 
h I=138 for N=81 read N=91 
- p=157,N=36 for I=80 read I =70 
326, p = 163,N = 143 for I =134 read I = 154 
p= 167,N=109 for I=36 read I = 35 
N=113 for I= 193 read I = 103 
N=161 for I=144 read I = 147 
I=147 for N=61 read N= 161 


ter I=162 fon N=21 read N=25 

327,p =173,N=21 for I=138 read I =38 
-he N=49 for I=61 read I=62 
- is N=57 for I=72 read I=92 
ter 328, p= 179, 1=79 for N=36 read N=33 
7, ’=181,N=16 for I=175 read I =172 
xX N=99 for I=192 read I = 102 
uld N=102 for I=79 read I=76 


329,p=191,N=91 for I[=9 read I[=99 
I=0 for N=0 read N=1 
I=172 fon N=18 read N = 138 
b= 193,N=42 for I=133 read I = 138 
330, p = 193, T= 115 for N=182 read N = 82 
331,p =199,N=11 for I=89 read I = 189 
N=5S6 for I=10 read I=20 
I=32 for N=12 read N=7 


jn- 
All the errata in the table of linear divisors of quadratic forms listed in 
ion MTAC, v. 1, p. 441 as appearing in the 1944 edition are present in this new 
of edition with the addition of the following new errata: 
342 
les _ pgp 
332 x? + 13¥ for 22+1 read 522 +1 
333 x? + 467" insert 167 
ll x? + 695 for 77 read 73 
944 xt + 71 for 27 read 237 
les, 
334 x? + 87¥ for 17 read 317 
x? + 8952 for 106 read 105 
335 x? + 939 insert 121 
x? + 959" for 36 read 363 
336 x? — 20y9 insert 83 
x? — 349 omit 117 
xt — 41¥ for 63 read 73 
x? — 437 for 7 read 97 
xt — 479° insert 65 


for 27 read 187 
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page 
336 x? — 51y* insert 175 

x? — 53¥ for 11 read 131 

x? — 55 insert 67 and 201 

x? — 58¥ for 67 read 65 
337 x? — 74¥° insert 253 
338 x? — 855% for 73 read 173 
339 x? — 101° for 378 read 373 

S. A. JOFFE 


515 West 110th St. 
New York 25 


UNPUBLISHED MATHEMATICAL TABLES 


EpitoriAL Note: The UMT Fine [see MTAC, v. 4, p. 101] now contains 
the following manuscript table: UMT 77[D].—E. C. Bower, Natural 
Circular Functions for decimals of a circle [MTAC, v. 3, p. 425]. For an un- 
published table concerning elliptic integrals see RMT 787. 


95[A].—INsTITUT FUR ANGEWANDTE MATHEMATIK, Eidg. Tech. Hochschule, 
Ziirich, Table of Binomial Coefficients. Manuscript in the possession of 
the Institute. 


This is a table of the exact values of the binomial coefficients, extending 
as far as 
1 
( "a = 10089 13445 45564 19333 48124 97256. 
H. RUTISHAUSER 
Ziirich, Switzerland 


96[A].—J. W. WRENcH, Jr., & L. B. Smitu, Values of the terms of the Gregory 
series for arccot 5 and arccot 239 to 1120 and 1150 decimal places, respec- 
tively. Mss. in possession of the authors. 


The table of individual terms of the Gregory series for arccot 5 gives, in 
the range 501D to 1150D inclusive, the first 820 terms of that series. Ex- 
clusive of zeros following terminating decimals, the total number of sig- 
nificant figures involved is 379,290. The companion table of terms of the 
series for arccot 239 consists of 1120D values of the first 235 terms. The 
total number of significant figures in this table is 131442. 

The sums of the positive and negative terms of each series are given to 
the corresponding degree of approximation. From these data approximations 
to arccot 5 and arccot 239 have been obtained correct to 1148D and 1119D, 
respectively, as confirmed by the ENIAC calculation of these numbers 
[MTAC, v. 4, pp. 11-15]. 

For the sake of chronological accuracy it should be mentioned that the 
final checking of the 1120D approximation to arccot 239 was completed by 
Mr. Smith on 24 July 1949, and the calculation of arccot 5 had been com- 
pleted by the writer the previous month except for checking the data beyond 
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850D. This checking had not been completed when the ENIAC computation 
was made in September, 1949, and subsequent comparison with these inde- 
pendently computed values of arccot 5 and arccot 239 revealed several 
discrepant figures which were found to be due entirely to seven errors of 
transcription of data and addition of terms in the previously unchecked 
portion of the calculation of arccot 5. Mr. Smith’s value of arccot 239 agreed 
perfectly to 1119D with the more extended approximation found by the 
ENIAC. By 6 October 1949 all discrepancies had been removed and the 
derived approximation to m agreed through 1118D with the value deter- 
mined by the ENIAC. 

As a by-product of the extension of arccot 5 from 850D to 1150D the 
earlier table of 2*, m = 1(2)1207 [MTAC, vol. 2, pp. 246, 374], has been 
extended to the range m = 1(2)1667. The same method of checking that was 
employed before, namely Fermat’s simple theorem, was retained. In addi- 
tion, the last entry in the table was multiplied by 2° to obtain 2”, which 
had previously been computed by H. S. UBLER but not published. The two 
values of this power agreed perfectly. 


J. W. Wrencnu, Jr. 
4711 Davenport St. N. W. 


Washington 16, D. C. 


97(C].—Rapio CoRPORATION OF AMERICA, Table of Logarithms to the Base 
2 for numbers 0.000 to 1.000, Mimeograph MS in possession of the RCA 
Laboratories Division, Princeton, N. J. Copy deposited in the UMT Fixe. 


The function log: p is tabulated to 4D for 
p = 0(.002).2(.01)1. 


The table was prepared to facilitate studies in the theory of information- 
handling. 


R. SERRELL 
RCA Laboratories 


Princeton, N. J. 


98[D, H].—Stmpney Jounston. Solutions of sinx = + cx. Manuscript in 
possession of the author, 81 Fountain St., Manchester 2, England. 


This manuscript gives a solution x of sin x = cx for c = 0(.002).78(.001) 
.927(.0005).9665(.00025)98575,9854(.0002).9892(.0001).9951(.00005).9979 
and of sinx = —cx for c = 0(.001).136(.0005).1815(.00025).201(.0002) 
.2054(.0001).2122(.00005).2148(.00002).217 and a large number of values of 
c in the immediate neighborhood of c = .217. There is also a table of solu- 
tions of sin x = (1 — c*)x for c = 0(.002).048. Throughout the entire table 
x is given to 8D with A, A’, A*, A‘. 


99[D, H].—SipNEY JOHNSTON. Solutions of x — sinx = c. Manuscript in 
possession of the author and NBSCL. 


This gives solutions to 8D of x — sinx = y*/6 for y = 0(.01)1.53. 
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100[F].—R. C. Doutuitt, Tables related to Euler's totient function. On punch 
cards in possession of the author, Dept. of Math., University of California, 
Berkeley. Tabulated manuscript in UMT FILE. 


We denote as usual by ¢(”) the number of numbers not exceeding m and 
prime to and define 
&(n) = z o(m) 
A(n) = [4 + 3n’*x] 
B(n) = ®(n) — A(n) 
C(n) = A(n) — ®(n — 1). 


The table gives all five functions form = 1(1)10000. The values of ¢() were 
taken from GLAISHER’s tables [MTAC, v. 1, p. 136], and the other functions 
were computed and checked by means of IBM 604 and 602A calculating 
punch machines. 

The functions B(m) and C(m), contrary to a conjecture by J. J. Syt- 
VESTER, are not always positive and are found to be negative or zero for a 
total of 20 times for B(m) and 18 times for C(m). [See also UMT 86[F], 
MTAC, v. 4, p. 29-30.] 


101(F].—D. H. Leumer, Tables of Ramanujan’s r(n). Tabulated manuscript 
and punched cards deposited in UMT FILe. 


The function r(m) is the coefficient of x"! in the expansion of the 24-th 
power of Euler’s infinite product 


(1 —x)(1 —x)(1 — x*)---. 
The tables give r(m) for m = 1(1)2500 as well as 
DX 7(n), x |r(n)|, LD {r(n)}? 
n<N n<N 


nsN 


for N = 10(10)2500. There is a separate table of r(p) for 1000 < p < 2500 
and a prime, giving also 6D values of |7(p) | p—"/*. The table was produced 
on an IBM 602A calculating punch. 


102[L].— CAMBRIDGE UNIVERSITY MATHEMATICAL LABORATORY. Table of 
1/I'(x + iy). Manuscript in the possession of the Laboratory. 


The table furnishes real and imaginary parts of the reciprocal of the 
Gamma function, 1/T'(x + ty), to 6D for x = — .5(.01).5 and y = 0(.01)1. 
This table was computed and printed by the EDSAC (MTAC, v. 4, p. 61-65) 
under the direction of J. P. STANLEY and involved nearly 20 hours of 
machine time. 
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103[L].—SipNEY JounsTon, Tables of Sievert’s Integral, manuscript in posses- 
sion of the author, photograph copy at NBSCL. 


This is a table of the function 


[ex (— A sec é)dt 


for A = 0(.5)10 and x = 0(x/180)x/2. Values are given mostly to 5S. Ex- 
plicit tabulation is not made beyond a value of x where the integral remains 
unchanged to 5S. 


104[V].—SrpneEy Kap tan, Tables of Velocity Functions Characterizing Flows 
Formed by Jets from Orifices. U. S. Naval Ordnance Laboratory Memo- 
randum, 87 p. Available only to government agencies and contractors. 


The basic mathematics governing the flow of incompressible fluids has 
been known for many years. Because a great amount of tedious computation 
is necessary, flows for only a few isolated cases have been calculated in the 
past. At the suggestion of G. BrrxHorr,' the Naval Ordnance Laboratory 
has calculated for the first time the flow patterns of an incompressible fluid 
from an orifice for four different angles of aperture: a = 0°, 15°, 45°, and 90°. 
In all, more than 2,000 points were calculated. 

The governing equations are 


(1) W = U+iV = In(g) — In (¢* — 2C¢ +1) 
(2) Z=xtiy=2' + iSe” 

(3) f= + in, Isi<1, O¢argge¢er 
(4) Z’= —-5°+C(W+ Inf) 





ll 


(5) z" in| { =p Ga 


jaar, a<argZ” Cr+ea 


where S = sin a, C = cosa, and values for a, U, V are as follows. 


a Range in U Range in V Number of Cases 
0 — 2.8(.2)2 0(x/20)x 525 

x/12 — 2.4(.2)3.2 0(x/20)x 609 

4/4 — 2.4(.2)1.6 0(x/20)x 441 

a/2 — 2.4(.1)1.6 O(x/20)x/2 451 


x and y are given to 4D and é and » are given to 5D. In each case the error is 
less than a half a unit in last place. 


1G, Birxnorr, & E. ZARANTANELLO, Jets, Wakes and Cavities, soon to be published. 
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Cannon, 225 Far West Building, National Bureau of Standards, Washington 25, D. C. 











AUTOMATIC COMPUTING MACHINERY 


TECHNICAL DEVELOPMENTS 


The Incorporation of Subroutines into a Com- 
plete Problem on the NBS Eastern 
Automatic Computer 


The construction of a modest-scale automatically-sequenced electronic 
digital computer at the National Bureau of Standards is nearing completion. 
This computer, called the NBS Eastern Automatic Computer (SEAC), is 
expected to be a useful tool both for numerical computation and for research 
in numerical analysis. In order to facilitate the exploitation of the SEAC, 
the Machine Development Laboratory of the Bureau has planned in advance 
key instruction routines for use on the machine. Groups of these instruction 
routines will be incorporated as subroutines in instruction programs govern- 
ing the solution of complete problems. It is the purpose of this paper to 
show how subroutines prepared in advance of problem solution, selected from 
a library of permanent subroutines, will be properly inserted by the com- 
puter in the instruction program relevant to a problem at hand. 

In order to follow this program it will be necessary to understand the 
principal features of the logical design of the SEAC. 

1. Memory.—The present routine is written to fit the initial model of 
the machine which has 64 acoustic lines, or tanks, each line storing eight 
“words.” A word consists of 48 binary digit positions, of which 45 are used 
to represent either a number or a four-address command, and the remaining 
three provide spacing between words. A number, N, of 44 binary digits 
followed by a sign digit, is stored as an absolute value with the proper sign 
attached. A plus sign is represented by “‘0’’ and a minus sign by “1.” The 
binary point of the number is located between the second and third positions 
from the left so that | N| < 4. 

In the operation representation, 10 binary positions are apportioned to 
each of four addresses, or memory locations. As there are only 512 memory 
locations in the initial model, nine binary digits are sufficient to specify any 
address (six digits to indicate the tank number and three digits to indicate 
the word within that tank). The first binary digit from the left of any address 
will always be zero and therefore will not be indicated here. The remaining 
nine binary digits will be represented by three octal digits. In addition, four 
binary digits represent, in coded form, the operations to be performed, 
designated herein by capital letters. The 45th digit from the left is again a 
sign indicator. 

The following important feature of the acoustic-line memory should be 
kept in mind: information sent to an address will replace the previous con- 
tent of that address. 

2. Command Code.—In Table I under the headings of each of the four 
addresses a, 8, y, and 6 the items which are located therein are listed. An 
address enclosed in parentheses indicates the content of the corresponding 
memory location; the address of the next command is abbreviated as ‘‘n.c.”’ 

The 45th binary digit from the left in a command is normally a zero 
(i.e., a plus sign); when a minus sign is coded in that position, the machine 
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Operation a 
Addition Augend 
Subtraction Minuend 
Multiplication Multiplicand 

(High-order) 

Multiplication Miultiplicand 

(Low-order) 

Multiplication Multiplicand 

(Rounded) 

Division Divisor 

Logie Extractee 
ransfer 

Algebraic ist Comparand 

i 
Absolute Value ist Comparand 

Comparison 
Input Order Any Number 
Input Order Any Number 
Output Order Any Number 
Output Order Any Number 


TABLE I 
8 7 
Addend Sum 
Subtrahend Difference 
Multiplier High-order 
oduct 
Multiplier Low-order 
Product 
Multiplier Rounded 
Product 
Dividend Quotient 
Extractor 
Number 
2nd Comparand n.c. or 
2nd Comparand n.c. or 
Odd Number Word from 
Input 
Even Number Word from 
Input 
Odd Number Word to be 
Recorded 
Even Number Word to be 
Recorded 


a 
s tar 
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nc. A 


n.c. 


n.c. 


nc. 


n.c. 


n.c. 


n.c. 


n.c. 


» oP 





Result of Operation 
() = @ + ®) 


(y) = @) - 


(vy) = the more significant half 
of (a)(6) 


(y) = the less significant half of 
(a)(6) 


(vy) = the more significant half 

of (a) @) increased by the first 

ty the less significant half 

the product 
(vy) = @)/(e), unrounded 
my poy to 
- are replaced 

responding digits of (a) 


(vy) = n.c., if (a) < (@) 
(6) = nc., if (a) > @) 


(y) = nc., if | (e) 

(6) = nc., # [Sls Bl 

o) =. = word from Input 

(y) to (y +7) = next 8 words 
from Input medium 

(y) = next word on Output 
medium 


(vy) to (y +7) = next 8 words 
on Output medium 


will stop automatically after executing the indicated operation. Commands 
can be manipulated in the arithmetic unit of the machine, since, to this unit, 
they are indistinguishable from numbers. 

3. Timing.—A minor cycle, the time needed for one word to pass a given 
point in the machine, is equal to 48 microseconds. The time consumed for a 
given operation, located at address ¢, can be obtained from the following 
formulae (where an underscored symbol represents the last octal digit of 


the indicated address) : 


1. For operations A, S, and L,¢ = (« — 8) + (8 — a) + (a —-y) + (vy — 9), 
where 8 minor cycles must be added to each difference which is alge- 
braically less than +1. 

2. For operations M, N, R, and D, the above sum must be increased by 40 
minor cycles if (a — y) > 5; otherwise, the sum must be increased by 


48 minor cycles. 


3. For operations C and K, two timings are possibie depending on the result 
of the comparison; the number of minor cycles is given by (¢-— 8) 
+ (8 — a) + (@ — 4), if (a) < (8), or (¢ — 8) + @ — a) + (a 9), if 
(a) > (8). In this case, 8 minor cycles must be added to the first two 
differences whenever their value is less than +1 and to the third differ- 
ence whenever its value is less,than +2. 
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The following table indicates the range of execution time for each command: 


Minimum Time Maximum Time 
Operation minor cycles microsec. minor cycles microsec. 
Cz oF 192 25 1200 
A, SL 4 192 32 1536 
M, N, R, D 48 2304 76 3648 


, About 2 seconds per word, at present 


4. Modification of Subroutines.—In 1947, SamueEL LUBKIN considered 
the problem of modifying subroutines which were to be used in a given 
problem to fit the available machine storage. He proposed the use of a Base 
Number Command to simplify the task of adapting any subroutine to its 
position in the memory. GOLDSTINE & von NEUMANN" have reported on the 
programming of such adaptation using the codes of the Institute for Ad- 
vanced Study machine. A similar routine prepared for use on the SEAC 
is presented herein. Because of the restricted storage capacity of this 
machine, an effort has been made to pack frequently-used routines in the 
internal memory as economically as possible. The present routine, in which 
the assumption is made that m subroutines are to be incorporated within the 
given program, occupies only 17 + cells. This number does not include, 
however, four temporary storage cells and the storage for nine constants 
which belong to a common pool (as explained in the next section). 

A general explanation of the method used in the modification of sub- 
routines for insertion in main instruction routines will perhaps be helpful to 
the reader. Each permanent subroutine is coded as though its first word were 
located at the address 0 100 000 000 (i.e., the binary equivalent of the octal 
address 400). All addresses in the subroutine which require modification, as 
introduced into the computer, are greater than 400. Thus, there will be a 
“1”” present in the second position of every address which must be modified 
and a zero in the corresponding position of all other addresses which are 
less than 400. 

In the present program, each of the required group of subroutines is 
inserted into the memory in the final locations it is to occupy. The modifying 
routine begins at position 050. All addresses of 400 or above contained 
within the subroutine are modified to fit the actual location of the sub- 
routine in the memory. In substance, the four digits corresponding to 
the second digit from the left of each address within such words of the 
subroutine as need modification are extracted. The resulting number is 
multiplied by d, which represents the difference between the address occu- 
pied by the first word of the subroutine under consideration and the number 
400 at which address this word was originally coded. The resulting product 
is added to the original word, and the addresses contained therein are thereby 
properly modified. 

Consider a specific example illustrating this procedure. As stated above 
the addresses will be written in octal form, although they are stored in 
binary form. Suppose the following subroutine operation, originally coded 
at address 400, is placed in the memory at address 070: 


Cell No. a B ¥ r Operation 
070 401 046 401 415 A 
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In order to modify (070), subtract 400 from 070 octally with the result, 
d = —310. Next, extract the second binary digit from the left of each 
address giving a “1” in the case of a, y, and 6 and a “0” in the case of £. 
Multiply each extracted digit by d = —310, giving —310 in the a, y, and 6 
positions of a memory cell and 000 in the 8 position, and add the result 
octally to (070). The modified operation is as follows: 


Cell No. a B Y é Operation 
070 071 046 071 105 A 


After all the subroutines in a program are modified, the memory cells 
occupied by the modifying routine are available for other uses. 

5. Conventions Used in Coding.—RBefore actually presenting the pro- 
gram, it will be necessary to explain some of the conventions used in the 
coding of this routine. An address enclosed within brackets is used to indi- 
cate the fact that the content of that storage location will vary. Braces are 
used in the ordinary algebraic sense. 

The first two memory tanks (i.e., memory locations 000 through 017) 
are retained for temporary storage in the modification routine; cell 007 
contains the instruction to be executed immediately after the routine is 
completed. Positions 020 through 047 serve as a common pool of frequently- 
used constants, which all subroutines employ. Of this pool, the present 
routine makes use of the following constants: 


Memory Cell Content 
020 + Zero 
023 400 (in binary form) occupying the a space, zeros occupying the 
remaining binary positions 
027 the number 2-*, representing a unit in the last position of a 
031 the number 2-# 
034 the number 2, representing a unit in the last position of y 
041 the number 27° 
042 2? — 2-%, representing 10 units.occupying the a space 
044 2-3 — 2-8, representing 10 units occupying the y space 
047 Zero (at the start of the routine only) 


The following notations will be used in the program: 


Notation Significance 
as; the address at which the j-th word of the i-th subroutine is located, 
where i = 1, 2, 3, - -nandj = 1, 2, 3, --- 
db; the number of words to be modified in the i-th subroutine 
Wij the content of cell a;; before modification 
Wij the content of cell a;; after modification 
Sk the binary digit of @,; in position k from the left, where k = 1, 2, 3, 
cae 
PROGRAM 
No. of 
Cell Oper- minor 
No. a B y 6 ation cycles Result of Operation 
050 070 044 O55 062 L 14 (055) = 010 013 a, 061 A 
062 070 031 O12 063 N_ S55 (012) = 2-*a,; (see footnote 2) 
063 041 O12 O12 067 D_ = 52 (012) = 2-*a;, 
067 012 042 #OS1 064 L 19 (051) = ag, 020 010 066 A 
064 012 O23 O12 O51 S ii! (012) = {ay — 400}2* =d 
051 [es] 020 010 066 A_ 12 (010) = ®j; 











AUTOMATIC COMPUTING MACHINERY 








No. of 

Cell Oper- minor 

No. a B ¥ & ation cycles Result of Operation 

066 010 065 047 054 L 10 (047) = 522° + 5422-19 + 5992-7? 

S3q2~*° = t 

054 047 012 013 05S M 63 (013) = td =c, correction 

055 010 013 [a;] 061 A 14 (a5;) = By +c = wy 

061 051 070 052/056 K_ 15/11 Test as to whether d; modifications have 
been made. 

052 027 O51 O51 060 A _ 10 If b; modifications have not been made, 
(051) and (055) are stepped up. 

060 034 O55 O55 O51 A 15 

056 [071] 020 070 057 A_ 16 All 5; modifications have been made; 
(070) = [{aa + 5; — 1}2-* + a2), 
where ¢ = 2, 3, ---, (a + 1). 

057 027 056 056 053 A_ 12 (056) is stepped up. 

053 020 070 050/007 K _ 19/20 Test as to whether » subroutines have 
been modified. 

Temporary and Constant Storage 

065 2° + 2-10 4 2-20 + 2-20 

070 ay + bh — 1}2-* + 42%; Lieu + bh - 1}2-% + a4;2-*] 

071 @a + be — 1}2-% + a2 

067 +2 {an +b, — 1}2-° + an2- 

070+ Zero 


In this program the commands are arranged in logical sequence. On the 
input medium, they would be arranged in numerical order of addresses of 
the memory cells containing them. The time consumed for each instruction 
to be modified is 6.5 milliseconds. 

MDL Staff 
NBSMDL 


1H. H. GoLpstTINnE, J. VON NEUMANN, Planning and Coding of Problems for an Electronic 
Computing Instrument, Part II, v. 3, Institute for Advanced Study, Princeton, N. J., 1948 
[MTAC, v. 3, p. 541-542]. 

2 Because the first two digits from the left in the result of low order multiplication on the 
SEAC are always zero, two instructions are required to put a, into the a position of 
memory cell 012. 


DISCUSSIONS 
A Note on “Is” and ‘‘ Might Be’ in Computers 


Recent press reports have aired the disagreement between the “brain” 
and the “antibrain”’ factions in the computer fraternity. The term “brain” 
is a bit fanciful and perhaps smells slightly like commercial advertising; 
however, I feel that some of the stories given to the press have been more 
definitely misleading. I would like to enter a plea for careful distinction 
between facts and fancies by scientific people who write and speak for the 
general public. 

A specialist in one field must be particularly careful in talking to the 
public about matters not in his own field of specialization. A careless or 
exaggerated remark made to experts in a field will do little harm because 
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the audience will make suitable allowance for “embroidery,” but a similar 
remark made to the general public may be taken seriously and may lead 
to real misunderstanding. 

At the present critical period in the development of science, when scien- 
tists are being forced somewhat reluctantly to take a hand in international 
affairs, there is much loose talk on the part of those who were formerly in 
undisputed control. Instead of excusing similar inaccuracies on the part of 
scientists, this seems to me to emphasize our responsibility to the public to 
present fair, accurate, and dependable statements. 

This note is frankly an appeal to those who issue frequent statements to 
the public press, to distinguish carefully between imagination and fact— 
between what “might be’’ and what “‘is.” I fear that these categories have 
been confused on several occasions. 

The analogy between computers and the brain or nervous system, for 
example, has been a useful aid to invention and, when properly qualified, 
to an understanding of the workings of a computer. Invention, as I imagine 
it, consists of selecting out of many random combinations of ideas some that 
serve a useful purpose. Any device, no matter how devoid of logic, that helps 
to form useful combinations should be accepted thankfully. In this role, by 
suggesting what might be, the “brain” analogy had been useful to me and 
probably to others, particularly in the formative period of computer develop- 
ment. Once a combination of ideas has been formed, however, it must stand 
on its own logical feet; one must then be careful not to distort facts to fit 
them into the analogy no matter how striking. 

According to an article in the Saturday Evening Post, Feb. 18, 1950, 
computers are subject to psychopathic states which the engineer in charge 
cures by a “shock treatment”’ consisting of the application of excessively 
large voltages. I would be interested to know whether this is a ‘“‘might be” 
or an “‘is.”” So far, I have been unable to find anyone who recommends quite 
so slap-dash a method of trouble-shooting. It would be too much like 
dropping a valuable watch on the floor to improve its time-keeping ability; 
granted that this treatment would work once in a while, it hardly seems the 
economical procedure in the long run. 

Another “might be”’ is the checking method that the reader of the public 
press would infer is generally used, with three parallel computing elements, 
all doing the same problem and voting on the correct answer. A majority 
vote decides the answer that the computer will print. There may be such a 
machine, despite the low economy of this checking method, but to date I 
have been unable to find this remarkably democratic device. 

The listing of ‘might be’s’’ could be extended to greater length. I should 
like to think, however, that I have touched upon a matter of concern to the 
entire scientific fraternity and that it is unnecessary to belabor my point. 
In all sincerity, I am urging more care and restraint in reporting scientific 
work to the public, particularly when the person reporting has only a second- 
ary interest in that work and is likely to be misled by incomplete information. 


G. R. Srrertz 
Burlington, Vermont 
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BIBLIOGRAPHY Z-XII 


1. E.G. ANDREwWs & H. W. Bone, “Use of relay digital computer,”’ Electrical 
Engineering, v. 69, 1950, p. 158-163. 


Use of two Bell Laboratories large-scale digital computers is described; 
discussion of process of setting up problem, decisions which must be made 
by attending mathematician, and those decisions which can be left to ma- 
chines are noted; samples of different types of problems which have been 
solved successfully by machines are given; and procedures used are discussed. 


Epitu T. Norris 
NBSMDL 


2. ANON., Description and Use of the ENIAC Converter Code, Ballistic 
Research Laboratories, Technical Note No. 141, Aberdeen Proving 
Ground, Maryland, November 1949, 23 pages. Mimeographed. 


In the first three years of its existence, the ENIAC, in spite of its fre- 
quent and often prolonged breakdowns, executed its tasks so speedily that 
the attending crew of programmers found it exceedingly hard to keep pace 
with its operation. The introduction of a set of coded orders by JOHN 
von NEUMANN and Aberdeen’s mathematical staff materially lessened 
the effort and time consumed by the attendants in supplying problems for 
the ENIAC. A new unit, the Converter, was incorporated which receives 
the code digits and initiates the corresponding operations. The machine 
operates on a one-address control system. 

The list of orders falls into essentially three categories. The storage orders 
call for the machine to perform the following operations: 


a) Transmit information from an IBM card to the eight Transmitter Groups, each 
holding a signed 10-decimal-digit number. 

b) Transmit information from the above eight groups and from an additional pair of 
similar nature to 20 Registers (18 of which are misnamed Accumulators, as pointed 
out by the authors themselves), each holding a signed 10-decimal-digit number. 

c) Transmit information from the three Function Tables to the above Registers. (A 
Function Table consists of 104 lines, each holding two signed six-digit numbers or 
a sequence of coded instructions. Each line bears a three-digit address.) 


The arithmetic orders call for the machine to perform the following oper- 
ations: 


a) Add a number to the contents of an Accumulator. 

b) Subtract a number from zero. 

c) Multiply two numbers with single precision, and add the product to a third number. 

d) Divide two numbers with single precision, with remainder available. 

e) Extract square root, with remainder available. 

f) Obtain the absolute value of a number. 

g) Obtain (1 — N)10*, where N is a number whose decimal point is k places from 
the left. 

h) Multiply a number by 10°, where — 5 < p < 5. (These shift orders are of two types, 
allowing for the retention of the shifted-out digits in another Register, or for their 
deletion. ) 


The control orders call for the machine to perform the following operations: 


a) Transfer Control—both conditionally and unconditionally. 
b) Count in order to facilitate iterations. 
c) Delay and stop. 
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The publication gives a careful explanation of each order but unfortu- 
nately exhibits some of the ills always accompanying the appending of an 
“afterthought.” For instance, the letters A and B are used to indicate the 
first two Transmitter Groups, as well as the two sides of each Function 
Table line. The letter X, often misprinted as x, represents the Multiplication 
Code as well as an arbitrary number within a register. A slightly more serious 
error is the wrong illustration of order 6L, at the top of page 10, which 
accompanies the correct verbal explanation. 

The interested reader would like to know whether there is any provision 
in the ENIAC for detecting overflow on addition and whether the multiplica- 
tion order gives a rounded product. This technical note was designed for use 
within the Ballistic Research Laboratories and was written primarily for 
issuing available information when the occasion demands speed. The exposi- 
tion of the codes is valuable in that it is sufficiently clear to enable any 
reader to attempt programming computations on the ENIAC. 


IDA RHODES 
NBSMDL 


3. ANON., Digital Computer Research at Birkbeck College, Office of Naval 
Research (London Branch), Technical Report OANAR-50-49, 12 De- 
cember 1949, 2 p. 


This report describes the design features and operating characteristics 
of the three digital computers being constructed under the direction of 
A. D. Boots of Birkbeck College, University of London. They are the auto- 
matic relay computer (ARC), the simple electronic computer (SEC), and 
the all-purpose electronic x-ray computer (APEXC). 

The ARC is a part relay, part electronic parallel machine using magnetic 
drum storage. Teletype tape is used as an input medium, and data is entered 
into the machine in binary form. A teleprinter or tape perforator will be 
used for output. Several simple problems have been run on this machine 
such as tabulation of prime numbers and factorization of large numbers. 
Although it is now temporarily dismantled, it will be reassembled and will 
use cross-bar type memory units. 

The SEC is a small, all-electronic, parallel calculator. It has the same 
type of magnetic drum memory as the ARC; teletype is used for input and 
output. 

The APEXC will be an evolution of the SEC using a magnetic drum 
store and employing about 800 vacuum tubes. The machine will store 1,024 
numbers of 32 binary digits. Input and output will be on magnetic tape. The 
machine will be used initially for making crystallographic computations. 


Epitu T. Norris 
NBSMDL 


4. ANon., The EDSAC Computing Machine, Cambridge University, Office of 
Naval Research (London Branch), Technical Report OANAR-43-49, 25 
November 1949, 5 p. 


This report describes the design features and present use of the EDSAC 
built by M. V. Wiikgs and his co-workers at the University Mathemati- 
cal Laboratory, Cambridge, England. [See MTAC, v. 4, p. 61-65.] In 
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order to gain experience with the programming of problems on the EDSAC, 
a number of simple computational problems have been run on the machine 
(e.g., the calculation of prime numbers, tabulation of the complex gamma 
function [UMT 102], computation of the number ‘‘e”’ to 200 decimal places, 
and calculation of the Arry integrals). 

Epita T. Norris 
NBSMDL 


5. WARREN S. McCut_ocu & JOHN PFEIFFER, “Of digital computers called 
brains,’ The Scientific Monthly, v. 69, 1949, p. 368-376. 


Comparison between known operations in the human brain and similar 
operations in large-scale computers. 


6. H. J. McSxrmin, “Theoretical analysis of the mercury delay line,” 
Acoustical Society of America, Jn., v. 20, 1948, p. 418-424, 4 figs., 
2 tables. 


The electromechanical aspects of problems presented by the use of 
mercury delay lines in computers are discussed. After a brief discourse on 
the theory of wave propagation for liquids, an analysis is made of the voltage 
developed by the piezoelectric pick-up crystal used in the delay line and of 
the distortion that might result when the carrier pulse is modulated. It is 
found that waves with slightly differing phase velocities can exist and may 
produce distortion effects. 

Epita T. Norris 
NBSMDL 


7. R. D. Ricotmyer & N. C. METROPpOLIs, ‘‘Modern computing,’’ Physics 

Today, v. 2, Oct. 1949, p. 8-15. 

Developments during the last decade of automatic computing methods 
with rapid progress involving engineering principles and techniques, mathe- 
matical methods, and logic of automatic computation are discussed. Also 
described are the functions of the ENIAC, methods of presenting problems 
to computing machines, possibilities of using machines for analytic pur- 
poses, and practical uses of computers in their present stage of development. 

EpitaH T. Norris 
NBSMDL 


8. A. E. Smita & C. V. L. Smitu, “Digital computers and their applica- 
tions,” American Society of Naval Engineers, Jn., v. 61, 1949, p. 137- 
168, diags. 


Presented herein is a general discussion of the design principles of 
digital computers and a description of existing machines and those under 
development. 


NEws 


Association for Computing Machinery.—Rutgers University was host to members of 
the Association at a conference on automatic computing machinery. The conference included 
general sessions on Tuesday morning and evening, March 28, plus two parallel groups of 
sectional meetings designated as Section A and Section B, March 28 and 29. 
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In Section A, the talks dealt with applications of computers, engineering descriptions 
of several complete machines, and engineering problems encountered in the design of 
machine components. Section B was concerned with the mathematical aspects of the design 
of automatic computers, the coding of complete problems to be solved on the machines, 
and mathematical descriptions of complete machines. 

In addition to the program of talks, the following exhibits were shown: the REAC by 
the Reeves Instrument Corporation, the RCA Linear Simultaneous Equation Solver by the 
RCA Laboratories, the IBM Card-Programmed Electronic Calculator by the International 
Business Machines Corporation, the small Binary-Octal Calculator by the Raytheon Manu- 
facturing Company, the MADDIDA (Magnetic Drum Digital Differential Analyzer) by 
Northrop Aircraft, Incorporated, the oscillograph in the computer field by Allen B. DuMont 
Company, and Zator card equipment by the Zator Company. 

The program was as follows: 


General Session, Tuesday, March 28, Dean Etmer C. Easton, College of Engineering, 
Rutgers Univ., presiding: 

Welcoming address by Mason W. Gross, Provost, Rutgers University. 

“High speed computing machines—a survey,” Perry Crawrorb, Research and De- 
velopment Board. 

“Machines of moderate cost,” GzorGE St1pitz, Burlington, Vt. 

Section A, Epwarp W. Cannon, National Bureau of Standards, presiding. 

“Applications of computing machines to the solution of management problems,” 
MARSHALL K. Woop, Department of the Air Force. 

“Engineering applications of electronic analog computers,” HERBERT ZAGOR, Reeves 
Instrument Corporation. 

“Applications of electronic computers to Census Bureau problems,” James L. Mc- 
PHERSON, FLORENCE K. Koons, RaLpo E. MULLENDORE, Bureau of Census. 

“Applications of the BINAC,” Joun W. Maucaty, Eckert-Mauchly Computer Cor- 
poration. 

Section B, Frep G. Fenper, Rutgers University, presiding. 

“The BINAC—a technical report,” James R. Weiner, Eckert-Mauchly Computer 
Corporation. 

“An analog series computer,”” Max G. SCHERBERG, Office of Air Research. 

“Preliminary design of the Mark IV,” Benyamin L. Moore, Harvard Univ. 

“Design of a low-cost computer” (read by C. V. L. Smita), Pau, L. Morton, Univ. 
of Calif. 

General Session, banquet, Houston Peterson, Rutgers Univ., toastmaster. 

“Automatic computing machinery of moderate cost,” Howarp H. Arxen, Director, 
Computation Laboratory, Harvard Univ. 

Section A, Wednesday, March 29, James J. Stave, Jr., Rutgers Univ., presiding. 

“A digital computer for solution of simultaneous linear equations,” SaMUEL LUBKIN, 
Electronic Computer Corporation. 

“The ANACON, a large-scale general-purpose analog computer,”” D. L. WHITEHEAD, 
Westinghouse Electric Corporation. 

“The IBM card-programmed electronic calculator,” CurHpert C. Hurp, International 
Business Machines Corporation. 

“The MADDIDA, general features,”” FLovp G. STEELE, Northrop Aircraft, Inc. 

Section B, Franz L. Att, National Bureau of Standards, presiding. 

“Planning and error consideration in the numerical integration of a difference equation,” 
Francis J. Murray, Columbia Univ. 

“Probability methods in the solution of elliptic partial differential equations,” Joun H. 
Curtiss, National Bureau of Standards. 

“‘A machine method for solution of systems of ordinary differential equations," RICHARD 
F. Curppincer & BERNARD DimspALE, Aberdeen Proving Ground. 








174 AUTOMATIC COMPUTING MACHINERY 


“The theory of digital handling of nonnumerical information and its implications to 
machine economics,’ CaLvin M. Mooers, Zator Company. 

Section A, Ernest G. ANDREws, Bell Telephone Laboratories, presiding. 

“New circuits installed in the Aiken relay calculator,” FREDERICK G. MILLER, U. S. 
Naval Proving Ground. 

“Digital computing machine components of universal application,’’ WILL1am S. ELLiotrt, 
Research Laboratories of Elliott Brothers, Ltd. (London). 

“Design features of a magnetic drum information storage system,” JoHn L. HILL, 
Engineering Research Associates, Inc. 

“The SB-256 electrostatic selective storage tube,’’ JAN A. RaycumMan, RCA Labora- 
tories. 

Section B, Eitis R. Ort, Rutgers Univ., presiding 

“Optical ray tracing,” DonaLp P. FEpER & BENJAMIN HANpDy, National Bureau of 
Standards. 

“Solution of matrix equations of high order by an automatic computer,’’ HERBERT F. 
MITCHELL, Jr., Eckert-Mauchly Computer Corporation. 

“'Theodolite reductions on the IBM relay calculators,” Mark LorxKin & C. E. JoHNson, 
Aberdeen Proving Ground. 

“The MADDIDA, design features,” DoNnALD E. ECKDAHL, Northrop Aircraft, Inc. 

The Institute of Radio Engineers.—At the 1950 National Convention held in New York 
City from March 6 through March 8, two sessions were devoted to high-speed computers. 

The first session on Wednesday morning, March 8, under the chairmanship of GEORGE 
R. Strsitz, dealt exclusively with digital computers. In a paper entitled, “Static magnetic 
pulse control and information storage,’’ AN WANG of Harvard University mentioned the 
possibility of developing a digital memory system using magnetic materials having rec- 
tangular hysteresis loops to control transfer of electrical pulses through the core by means 
of the property of residual magnetism. The incorporation of unidirectional current devices 
permits a static, high-efficiency pulse power distribution system. RaLpa SLutz of the 
National Bureau of Standards next discussed the development of a generalized procedure of 
designing diode circuits for pulse gating, with particular attention to the extremely high 
reliability required for electronic digital computers. In the third paper, ‘‘Marginal checking 
as an aid to computer reliability,” NorMAN H. Taytor of M.I.T. pointed out the source of 
error in digital computing and pulse communication caused by deterioration of machine 
components. Marginal checking varies voltages in logical circuit groups, inducing inferior 
parts to cause failure, while a test program or pulse transmission detects and localizes poten- 
tial failure. By automatically carrying out this marginal checking, a digital computer can 
act as its own detector. ‘The development of the California Digital Computer” by Davip 
R. Brown and Paut L. Morton, University of California, Berkeley, presented a description 
of a low-cost general-purpose electronic digital computer having a compact magnetic-drum 
memory with a capacity for 10,000 ten-decimal-digit numbers. This large memory capacity 
permits much computation without the necessity of automatic access to input. The average 
access time is eight milliseconds. Machine input is accomplished on punched paper tape 
read photoelectrically. Binary-coded decimal numbers pass serially through the arithmetic 
unit on four parallel channels, and single-address operations (including division and square 
rooting) will be performed at the rate of about 60 per second. In the final talk of this session 
JoseF Katz, University of Toronto, mentioned a new class of switching tubes designed to 
reduce the complexity of electronic switching circuitry. The tube electrodes (one for each 
input and output channel) are designed so that particular combinations of input voltages 
result in current flow to the corresponding combination of output electrodes. The total 
current, and hence the cross section of the electron beam, is restricted by a limiting resistor 
used to focus the beam. 

The second computer session on Wednesday afternoon, March 8, under the chairman- 
ship of Witt1am H. Hucerns, was concerned with information analysis and computing. 
The REAC (Reeves Electronic Analog Computer) was described by H. I. ZaGor, Reeves 
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Instrument Corporation, New York, as a flexible and practical tool for solving ordinary 
linear and nonlinear differential equations. In addition, the speaker illustrated the capa- 
bilities of the machine components and showed the various techniques involved in solving 
such representative problems as flutter, electron flow, automatic pilot design, Fourier 
analysis, engine control, integral and boundary value equations. ‘An electronic storage 
system” by E. W. Bivans and J. V. Harrincton, Air Force Research Laboratories, Cam- 
bridge, Mass., described the RCA Radechon (a barrier grid storage tube) which is being 
used in a digital storage system. The secondary collector system is not used; instead, the 
reading signals are measured at the back plate, the same electrode on which the write signals 
are impressed. Deflection voltages are generated by a weighted addition of the plate voltages 
of a binary counter. Read and write operations are asynchronous with a 12-y sec minimum 
time between operations. An application of the theory of correlation functions used in the 
determination of the transfer functions of linear and nonlinear systems was presented next 
by J. B. Wiesner and Y. W. Lee, M.I.T. In the paper, ‘Measurement and analysis of 
noise in a fire-control radar” by R. H. E1sENGREIN, Sunstrand Machine Tool Company, 
Rockford, Illinois, an ‘instantaneous subtraction” method of optically measuring radar 
noise in order to analyze the unwanted portion of the signal return from an airborne target 
was discussed. Finally, H. E. SiIncLteton, M.I.T., described a new electronic correlator which 
is capable of accepting inputs covering a wide frequency range and which evaluates correla- 
tion functions for arguments from 0 to 0.1 seconds. In order to obtain a high degree of 
accuracy and stability, the signals are sampled and converted to binary numbers, and the 
storage and computation are carried out digitally. 


OTHER AIDS TO COMPUTATION 
BIBLIOGRAPHY Z-XII 


9. D. P. Apams & H. T. Evans, “Developments in the useful circular 
nomogram,” Rev. Sci. Instruments, v. 20, 1949, p. 150-154. 


A circular nomogram for W = UV with the U, V scales on the circum- 
ference and the W scale on a diameter is described, with a discussion of the 
most advantageous choice of scales. 


10. J. A. Bronzo & H. G. Couen, “‘Note on analog computer design,” Rev. 
Sci. Instruments, v. 20, 1949, p. 101-103. 


This note proposes that partial differential equation problems be at- 
tacked by transforming them into difference-differential systems and then 
into systems of linear equations. The characteristic roots of the matrices of 
the latter are to be investigated by a change of coordinate system. The 
choice of the new coordinate system is not specified but is dependent on the 
ingenuity of the investigator. 


F. J. M. 


11. J. H. Fevxer, “Calculator and chart for feedback problems,” I.R.E. 
Proc., v. 37, 1949, p. 1204-1206. 


The chart consists essentially of constant-magnitude and constant-phase 
loci of z in the complex plane of y = z/(1 + z). A calculator based on this 
chart is available commercially. 

Lotri A. ZADEH 
Department of Electrical Engineering 
Columbia University 
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12. Pattie G. HuBBARD, “Application of electrical analogy in fluid me- 
chanics research,’’ Rev. Sci. Instruments, v. 20, 1949, p. 802-806. 


Fluid flow problems are solved approximately by the measurements on 
the electric potential in a suitably shaped electrolytic tank. The electric 
potential is analogous to the velocity potential of the flow, inasmuch as both 
satisfy LAPLACE’s equation and analogous boundary conditions can be ob- 
tained. Two techniques are described. In one of these, lucite models are 
immersed in a tank to simulate the flow pattern around the model and in the 
second, a tank segment with triangular cross sections was used to determine 
boundary pressures for a flow confined by various surfaces of revolution. A 
graphical comparison is made with experiment and the discrepancies 
analyzed. 

F. J. M. 


13. J. Lyman & P. V. Marcuett1, “A device for facilitating the computa- 
tion of the first four moments about the mean,” Psychometrika, v. 15, 
1950, p. 49-55. 


The device consists of a box with a number of endless tapes in it, ar- 
ranged so that a chosen line on each tape can be read through a window in 
the box. A given distribution is divided into equally spaced intervals. Each 
tape corresponds to an interval x and each line corresponds to a possible 
frequency f for this interval. On this line f, fx, fx’, fx*, and fx* are typed and 
' thus when each tape has been adjusted to the correct frequency for the given 
distribution, one can immediately read off the various terms which appear 
in expression for the first four moments. 

F. J. M. 


14. M. W. Maxowsk1, ‘Slide rule for radiation calculations,’’ Rev. Sci. 
Instruments, v. 20, 1949, p. 876-884. 


This paper describes a special slide rule for making calculations connected 
with PLANCK’s radiation formula. By setting the cursor to a given tempera- 
ture there may be found on the stock the total radiant energy and photon 
flux densities and the wave length of their maxima per unit wave length and 
per unit frequency as well as the values of their maxima. By setting the slide 
for a given temperature and the cursor for a given wave length, values of 
the ratios of the densities at any wave length to their maxima and of the 
integrated energy or number of quanta below or above that wave length 
may be read on the slide—WIENn’s law makes this possible. It is then simple 
to obtain actual densities and total fluxes within wave length ranges. Scales 
relating wave length, wave number, and electron volts are on the back of 
the stock. The accuracy is estimated by the inventor at from a few tenths 
of a per cent to five per cent, being of the order of one per cent in most ranges. 

L. H. THomMas 
Watson Scientific Computing Laboratory 
612 116th Street 
New York 27, N. Y. 
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15. R. E. Meyverott & G. Brett, “Small differential analyzer with ball 


catriage integrators and selsyn coupling,” Rev. Sci. Instruments, v. 20, 
1949, p. 874-876. 


This article describes a small differential analyzer built at Yale Univer- 
sity around three ball carriage type mechanical integrators. The advantages 
of this type of integrator are that it has a good output torque, the input and 
output shafts are mechanically convenient, and the units are readily avail- 
able as war surplus equipment. Coupling between integrators and input or 
output devices is by means of selsyns instead of shaft couplings, which 
makes the equipment very flexible and minimizes the set up time for new 
problems. The selsyns rotate at high speed (2500 to 1 to the integrator) to 
eliminate selsyn error. The analyzer also includes three input or output 
devices (called ‘function units’) and an “‘adder” consisting of two selsyns 
coupled to a mechanical differential. 

A test problem set in to generate a pure sine wave gave an increase of 
amplitude of 0.2 per cent per cycle, which compares favorably with results 
obtained with more elaborate machines. The device has proved satisfactory 
and convenient for the solution of problems with a nominal accuracy of one 
per cent. 

SipNEY GODET 
Reeves Instrument Corporation 
New York 28, N. Y. 


16. Cyrit STANLEY SMITH, “‘Blowing bubbles for a dynamic crystal model,” 
Jn. Applied Physics, v. 20, 1949, p. 631. 


This communication deals with the development of a certain type of 
analogue solution for problems in crystal structure and describes an im- 
proved method, using a nozzle, for obtaining the bubble rafts used by BRaGG 
and Nye for metal crystal study. (Cf. R. Soc. London, Proc., v. 190 A, 
1947, p. 474-481.) 


F. J. M. 


NOTES 


116. ON FINDING THE SQUARE ROOT OF A COMPLEX NUMBER.—There are 
two classical methods used for finding x and y in 


(1) (a + bi)! = x + yi. 
The first of these is based on DEMorvre’s theorem. One computes in 
succession 

R = (a + 8°)! ¢@ = arctan (b/a) 

x = Roos $(¢ + 2Kz) y = Rsin }(¢@ + 2Kzn), (K = 0,1). 

The second method is derived from purely algebraic considerations. 
From (1) it follows that 
(2) 2—y =a, 2xy = b 

R=x#+y=(¢+ 5)! 


x = + ((R + a)/2)! 
(3) y = + ((R — a)/2)'. 
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The signs are determined in accordance with (2). The second method avoids 
the use of trigonometric lookups. In addition, the first method is reducible 
to the second. Hence, the second method is the better one to use. However, 
it, too, has some serious shortcomings. When ) is small, either x or y suffers 
from a loss in significance. How does one prevent such a loss? 

GARRETT BrrKHOFF has suggested a third procedure. Define 


(4) t= }(R+ |a)). 

Then x and y are determined as follows: 
If a>d0, x = #, y = $br-3 
If a <0, x = $bt-}, y = #. 
If a= 0, x = y = (b/2)}. 


Of course in all cases, the signs of x and y are determined in accordance 
with (2). 

The numerical analyst should prefer the third method for two reasons: 
(a) accuracy is preserved, (b) the square root in (3) is replaced by a division, 
a great time-saving advantage. 

Example: Find the square root of —7 + .03%. 

The second method gives only 


x = .00566 95 y = 2.64575 7386 
while the third method gives 
x = .00566 94540 7 y = 2.64575 7386. 


Should any of the readers find an exposition of this technique, I should 
be glad to know about it. Birkhoff and I believe it to be new. 
SIDNEY KAPLAN 
Bureau of Employment Security 
Washington 25, D. C. 


117. TABLE FoR So_vinG CuBic EquaTions.—Judging by the title of 
the following 75-page pamphlet published more than sixty years ago one 
would never suspect that in it might be any table of mathematical interest: 
JoHN BorDEN, Some Higher Plane Curves, Chicago, E. T. Decker, 1888. 
Curves of the third, fourth, and fifth degrees are considered; but on pages 
55-75 is a table of the roots of m* — 3m + C = 0 for values of C between 
+2 and —2. Each column headed respectively 1, 2, 3, contains the values 
of a root in numerical order. The value of C next to the right of the column 
of roots has its proper value for C directly opposite and to the right of it. 
In the first column headed 3, the values are successively V3, 1.733(.001)2, 
and corresponding to each of these values is a 9D value for C. In the column 
headed 1, the values range 0(.001)1; in the column headed 2, the range is 
V3, 1.731( —.001)1. The author remarks: ‘‘The value of each root has its 
proper value for C directly opposite and to the right of it.” ‘When the value 
of C is positive, roots nos. 1 and 2 are positive and no. 3 is negative: and 
when C is negative nos. 1 and 2 are negative and no. 3 is positive.’’ Corre- 
sponding to the root 3 = 1.949, there are the corresponding root 1 = .945 
and the root 2 = 1.055, followed by 54 other pairs .946(.001).999 and 
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1.054( —.001)1.001. To each of these 109 roots is given the corresponding 
value of C to 9D. Then comes the final entry in the table: root 3 = 2, C = 2; 
root 1 = 1, C = 2; root 2 = 1, C = 2. 

» & A 


118. NOTE ON THE PAPER “On a definite integral” By R. H. Ritcuie 
[M7AC, v. 4, p. 75-77].—The value given by Ritchie for the integral 


f(x) = f: (u + x) exp (— u*)du 


is not free of ambiguity, since the definition of Ei(y) given in the paper does 
not make sense for y = x* > 0. The difficulty arises over the inversion of the 
Laplace transform (L.T.) p-' In (p — 1). 

Moreover, f(x) can be evaluated by a routine process, since it is simply' 
the repeated L.T. of exp (— u*). The L.T. of exp ( — w’) is? 


3x1 exp (47) erfc (44) 
where 


erfc zs = 2x73 f exp (— u*)du. 


Next, the L.T. of this can be obtained by means of CampBELL & Foster,” 
formula 959.5. The result is incorrectly rendered in Cossar & Erpéty1;* the 
corrected form is 


— 4exp (— x)Ei(— x*e*) + fix exp ( — x”) erfc (xe**!*). 
This can be expressed in terms of real valued functions if we use the function‘ 


Ei(s) = 4Ei(— ze*) + $Ei( — 2e-**). 


The result is 


f(x) = exp (— x*) [= f ; exp (u*)du — IEi(e) |: 


1D. V. Wipper, The sogioes Transform, Princeton, 1941, Chapter VIII, 1. 
2G. A. CAMPBELL & R. Foster, Fourier Integrals for Practical Applications, New 
York, 1948, no. 903.0. 
. Cossar & A. Erv& yi, Dictionary of Laplace Transforms, London, 1944-1946, p. VI, 
76. ie entry 3, read (2aix4)- ‘instead of ¢ ixt(2a)—; in entry 4 read erfc [p/(2a)] in "a of 
erfc (pat), and =~ Ei in place of Ei. 
*E, JAHNKE & F. Empe, Tables of Functions, Leipzig, 1938, p. 2. 


119. Harry BATEMAN BiBLioGRaAPHY.—In MTAC, v. 3, p. 141-142 we 
presented a complete summary of published biographical sketches of Bate- 
man, and lists of his publications, as well as a supplement to such lists of 
material published in the Educational Times (E.T.) and Educational Times 
Reprint (E.T.R.). Miss C. BrupNo, in the Applied Mathematics Branch, 
Mechanics Division, of the Naval Research Laboratory, Washington, has 
noted the following five additions to this latter list of problems proposed: 
14943, E.T., v. 54, 1901, p. 328; sols. by Bateman and another, £.T.R., v. 1, 
1902, p. 98-100. 14975, E.T., v. 54, 1901, p. 423; sols. E.7.R., n.s. , v. 2, 1902, 
p. 111, v. 3, 1903, p. 29. 15119, E.T., v. 55, 1902, p. 233; so 1. v. 55, 1902, 
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p. 516, also E.7.R., n.s., v. 3, 1903, p. 110-111. 15221, E.T., v. 55, 1902, p. 
438; sol. E.T.R., n.s., v. 4, 1903, p. 88. 16009, E.T., v. 59, 1906, p. 270; sol. 
E.T.R., n.s., v. 11, 1907, p. 57-61. 

R. C. A. 


QUERY 


35. W. THIELE’s TABLE.—According to KaysER’s VOLLSTANDIGES 
BUCHERLEXICON and Hinricn’s Katalog, Tafel der Wolfram’schen hyper- 
bolishen 48 stelligen Logarithmen. Bearbeitet und erweitert von W. Thiele, was 
published in two printings, by different companies, 118 p., at Dessau: (1) 
1905; (2) 1908. This publication appears to be extremely rare. There are 
copies of (2) in libraries of Harvard University, of Mr. C. R. CosEns, 
Cambridge, England, of Brown University (film), and in the John Crerar 
Library, Chicago. It was Cosens who, in 1939, directed my attention to (1). 
In what libraries is a copy of (1) located? 

R. C. A. 


QUERIES—REPLIES 


45. THe INTEGRAL Jf,’ e~4*dé (Q 19, v. 2, p. 196).—Attention is called 
to a manuscript table of this integral in UMT 103. 


46. Pitiscus TABLEs (Q 29, v. 3, p. 398; QR 40, p. 498-499, 42, p. 562- 
563).—For nearly 40 years I’ve had in my library BAsm_ ANDERSON & R. T. 
RICHARDSON, Catalogue of the Books and Tracts on Pure Mathematics in the 
Central Library, Newcastle-upon-Tyne, Newcastle-upon-Tyne, 1901. This 
library contains many old valuable books, but I did not earlier think to 
check for its possible Pitiscus items. On p. 33 it is indicated that the library 
owns copies of both 19966a (the Canon) and 19967 (the English Trigo- 
nometry, 1614), which therefore supplements the information we assembled 
[MTAC, v. 3, p. 499]. The library has also a copy of the 1612 German 
edition of the Trigonometry 6 [v. 3, p. 391]. 

R. C. A, 
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p. 54, 1. 3, for 4-17, read 3-16. 
p. 406, 1. 10 and 11, interchange cos }kx(k =1, 2, 3) and cos }(2k+1)x(k=0, 1, 2, 3). 
p. 21, 1. 1 and 2, for Asymptotic distribution of range from that of reduced range 
e distribution of the range. 1. 7 for 193-196 read 395-396. 1. 9 delete .1(.1).9(.01). 1. 
40D read ~. 
. 23, 1. —19, for 43-126 read 113-126. 
. 29, 1. —13, for p. xx read p. 11-15. 1. —12, for Jardan read Jarden. 
= 5 -1, for G10; /an read aj Qi: (a1;/an). 
8, 1. 25, for 1938 read 1935. 
2, 1. 14, for 548-553 read 948-953. 
1. 
1. 
1. 
1. 
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V.4 
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5, for SD read 4D. 

7, for U read u. 

-23, for 37-41 read 33-41. 
-23, for 332-371 read 322-371. 








